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ABSTRACT 

The  application  of  LORAN-C  in  the  hyperbolic  mode  as  a 
positioning  system  for  hydrcgraphic  surveys  was  investi- 
gated. Observed  LORAN-C  time  differences  from  a  field  test 
conducted  in  Monterey  Bay,  California  were  compared  to 
calculated  time  differences  determined  from  geographic  posi- 
tions based  on  a  microwave  positioning  system.  Four  methods 
were  used  to  determine  th?  calculated  time  differences.  The 
first  three  methods  were  (1)  applying  only  che  seawater 
Secondary  Factor,  (2)  computing  the  time  difference  based  on 
a  Semi-Empirical  TD  Grid,  and  (3)  applying  ASF  Correctors 
from  the  DMAHTC  LORAN-C  Correction  Table.  The  final  method 
applied  multiple  observed  ASF  Correctors  at  five  minute 
latitude      and   longitude      intervals.  By   applying      multiple 

observed  ASF  Correctors,  which  was  the  most  accurate  method, 
a  38.3  me-cer  1  drms  with  a  lane  offset  of  3  to  12  me-ers 
using  the  9940  X-Y  LORAN-C  combination  was  obtained.  Based 
upon  rhe  results  presented,  it  may  be  possible  to  use 
LORAN-C   for   hydrographic   surveys   at    scales   of    1:80,000. 
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I.       INTRODUCriON 

A.       USE    OF    LORAN-C 

In  recent  years  there  has  been  an  increasing  usage  of  a 
LORAN-C  receiver  and  the  LORAN-C  network  as  rhe  primary 
horizontal  control  for  such  scientific  studies  as  deep  ocean 
dumpsites,  marine  fisheries  studies  [Rulon,  1979],  bathyme- 
tric  surveys,  and  recently,  a  recDnn aissance  hydrographic 
survey.  Examples  of  bathymetric  surveys  conducted  by  the 
National  Oceanic  and  Atmospheric  Administration  (NOAA)  , 
which   have    used   Loran-C   for    positioning    are: 

1)  Su-100-1-79      Gulf    of    Alaska   [NOAA    H-9822,     1979], 

2)  SU-100-2-79      Gulf    of    Alaska   [NOAA    H-9823,     1979], 

3)  S-D902-WH-82    U.S.    West   Coast    [NOAA    Ship   Surveyor, 
1982]. 

Recently  an  attempt  was  made  ♦:o  use  LORAN-C  as  the 
sounding  position  control  for  a  reconnaissance  hydrographic 
survey  S-K902-WH-82  [NOAA,  1982].  This  was  a  special  survey 
conducted  by  the  NOAA  Ship  Whiting  in  May  1982.  Special 
surveys  are  field  examinations  of  very  limited  extent  or 
scope  and  frequently  require  unique  survey  or  data  collec- 
tion procedures  [Umbach,  1976].  The  purpose  of  this  project 
was  to  verify  the  existence  and  extent  of  reported  shoaling 
in   three   safety    fairways    in    the    Gulf   of    Mexico   [NOAA,    1982]. 

The  use  of  Loran-C  as  a  positioning  system  for  basic 
hydrographic  surveys  has  been  very  limited  due  to  the  abso- 
lute accuracy  of  the  long  range  system.  A  basic  hydro- 
graphic  survey  is  defined  as  a  survey  which  is  so  complete 
that  it  need  not.  be  supplemented  by  other  surveys.  "It  must 
be  adequate   to   supersede    for   charting    purposes   all    prior 


surveys"  [Umbach,  1976].  Variables  which  affact  the  accu- 
racy of  LORAN-C  are  signal  propagation  variations,  weather, 
and  sky  waves.  The  affects  of  weather  and  sky  waves  on 
LORAN-C  propagation  are  best  described  by  Samaddar  [1980] 
and  the  American  Practical  Navigator,  [DMA,  1977]  respec- 
tively. 

Signal  propagation  variations  are  due  to  the  phase 
retardation  of  the  signal  as  it  passes  over  an  all  seawater 
path,  over  land  paths,  or  partial  seawater-land  paths  as 
compared  to  free  space.  Table  I  sammarizes  phase  retarda- 
tion   changes   [Mortimer,         1978].         Errors    due   to      an    all    sea 

TABLE    I 
Phase   Retardation    or  Lag   of    Radio    Waves 


Propagation    Path 


ReDr=s  entation 
Propa  gation 
V9l ocity 
(kffl/secf 


Vacuum 


299792.5 


Direct    wave   through        299  691 
earth's   atmosphere 


Ground   wave  over 

sea    water 

Ground   wave   over 
rugged   mountains 


299560 


298  899 


Difference    in    Phase 
Lag   at    500    km 
Compared    with 
Wave   in    Line 
Abcv5     (m) 


170 


220 


1,300 


water  path  are  known  as  the  Secondary  Factor  (SF)  and  errors 
due  to  a  land  path  or  mixsd  path  are  known  as  the  Additional 
Secondary  Factor  (ASF)  [Speight,  1982]. 
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ASF  Corrections  in  ths  LORAN-C  system  can  b9  as  large  as 
plus  or  minus  four  microseconds,  which  is  600  meters  on  t.he 
baseline.  In  other  areas  with  the  same  LORAN-C  coverage, 
these  corrections  may  be  much  larger  due  to  the  expansion  of 
the  distance  between  adjacent  hyperbolic  lines  of  position. 
For  example,  at  32°  N  and  80°  W,  using  lattice  pair  9960-X, 
a  four  microsecond  (jisec)  error  will  offset  the  9960-X  line 
of   position   approximately   2438   meters    [Speight,    1982]. 

B.       APPLICATION    OF    ASF    CORRECTORS    FOR    NAVIGATION 

To  compensate  for  the  Loran-C  positional  errors  caused 
by  the  ASF  Correctors,  the  Secretary  of  Transportation 
tasked  the  Defense  Mapping  Agency  aydrographic/Topographic 
Center  (DMAHTC) ,  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  ,  and  the  United  States  Coast  Guard 
(USCG)  with  the  job  of  determining  and  applying  the 
Additional  Secondary  Correctors  for  each  Loran-C  chain. 
This  task  was  published  in  the  Department  of  Transportation 
(DOT)  National  Plan  for  Navigation  in  the  July  19,  1974 
Federal  Register.  These  corrections  should  provide  95^ 
assurance  that  a  vessel  could  fix  its  position  to  a  pred- 
icted accuracy  of  1/4  nautical  mile  (NM)  within  the  U.  S. 
Coastal  Confluence  Zone  (CCZ)  and  the  Great  Lakes.  The  CCZ 
is   defined    as: 

"the  seaward  approaches  to  land,  the  inner  boundary  of 
which  IS  the  narbor  entrance  and  the  cuter  boundary  of 
which  is  50  nautical  miles  offshore  or  the  edge  of  the 
Continental  Shelf  (100  fathom  contour)  whichever  is 
greater. " 

The      1/4-NM        accuracy      reguirement      also        affects      the 

nautical  chart.  The   National      Ocean    Survey       (NOS) ,      which 

publishes    charts   for   the    CCZ,      engaged    in    a    program    with    the 

USCG    and   DMAHTC    to    provide    the   coastal    navigator    with   charts 

overprinted   with   lattices   which      meet    1/4-nm    accuracy.         The 

OSCG,      as      operator   of      the    LORAN-C      radionavigation    system. 
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conducts  surveys  to  ensure  that  LORAN-C  coverage  exists 
within  the  CCZ  and  will  be  reponsible  for  the  verification 
of    1/U-NM      accuracy    for      all   coastal      LORAN-C   service.  In 

conjunction  with  NOS ,  it  assists  in  surveys  of  coastal 
waters  of  the  United  States  to  allow  production  of  LORAN-C 
charts  based  on  observed  field  data  to  meet  the  standards 
set    forth  above   [Speight,    1982]. 

DMAHTC,  for  LORAN-C  civil  need,  prepares  grid  predic- 
tions from  its  data  base.  Based  on  analysis  and  verifica- 
xion  of  the  predicted  grid  from  a  USCG  and/or  NOS  survey,  it 
produces  revisions  to  the  initial  grid  predictions  [ Speight, 
1982].  At  present,  DMAHTC  has  provided  NOS  with  ASF 
Corrected  LORAN-C  Lattices  which  are  overprinted  on  the  NOS 
Charts.  Each  chart  with  ASF  Correctors  applied  contains  one 
of  the   following   notes: 

"The  LORAN-C  lines  of  oosition  overprinted  en  this  chart 
have  been  prepared  for  use  with  groudwave  sianals  and 
are  presently  compensatsd  only  for  t heoretical' prooaga- 
tion  delays,  wnich  have  not  yst  been  verified  by 
observed  data.  Marinsrs  are  cautioned  not  to  rely 
entirely  on  the  lattices  in  inshore  waters.  Skywave 
corrections  are    not    provided". 

or 

"The  LOEAN-C  lines  of  position  overprinted  on  this  chart 
have  been  prepared  for  use  with  groundwave  signals  and 
are  compensated  with  propagation  delays  cofflDu.ed  from 
observed  data.  Mariners  are  cautioned  not  to  rely 
entirely  on  the  lattices  in  inshore  waters.  Skywave 
corrections  are   not    provided"    [Speight,    1982]. 

Presently,         all      of      the      NOS      Charts      of      1:80,000      to 

1:120,000    scale   covering      the  east    coast.      Gulf      coast,      and 

Great    Lakes   show   LORAN-C   lattices   that    have   been    compensated 

for    Additional   Secondary   Factors.         Most      of   'h=    lattices    on 

these    charts      have   been  constructed      from    DMAHTC      data   tapes 

that      provide   adjusted      LORAN-C   readings      for      each    rate      at 

every      five      minutes      of   latitude      and      longitude.  A      few 

lattices   were   constructed   using    a      single    ASF   Correction    for 

the    entire      chart    area.  Five    minute      data   tapes      were    not 
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furnished  by  DMAHTC  for  constructing  lines  of  position  for 
LORAN-C  rates  on  the  West  Coast  Charts.  On  these  charts  a 
single  average  ASF  Correction  was  used  to  adjust  each 
lattice   [NOAA,    Marine    Chart    Division,    1982]. 

In  addition  to  supplying  corrected  LORAN-C  lattices  for 
nautical  charts,  DMAHTC  prepares,  distributes,  and  periodi- 
cally updates  unclassified  ASF  LORAN-C  Correction  Tables 
[Speight,  1982].  The  ASF  Correction  Tables  are  for  preci- 
sion navigation,  utilizing  digital  computers  to  convert 
LORAN-C  time  differences  to  geographic  coordinates  [Speight, 
198  2].  Presently,  the  ASF  correctors  found  in  the  LOEAN-C 
Correction  Tables  were  determined  using  theoretical  propaga- 
tion delays.  ASF  correctors  listed  in  the  tables  are  going 
to  be  updated  with  observed  data  and  reprinted  the  first 
quarter    of    1983    [Wallace,    1982]. 

C.       APPLICATION    OF    ASF   TO    HYDROGRAPHIC    POSITIONING 

Schr.ebele  [1979]  investigated  the  possibility  of  using 
Loran-C  as  an  electronic  positioning  system  for  hydrographic 
surveying.  He  concluded  that  in  Monterey  Bay,  California  a 
single  Additional  Secondary  Factor  (ASF)  applied  to  offshore 
lines  of  position  gave  a  root  mean  sguare  error  (drms)  of  66 
meters  for  the  West  Coast  99U0  Y-W  pair  and  a  predicted  U2 
meter   drms    error   for    9943   X-Y  rates. 

The  42  meter  predicted  drms  is  larger  than  Nelson's 
[General  Electric  Co.,  1979]  findings  in  San  Francisco  Bay. 
He  demonstrated,  in  a  dynamic  mode,  that  the  precision  of 
LORAN-C  was  60.8  meters  2  drms  (30.4  meter  1  drms)  with  a 
worst  case  of  71.2  meters  2  drms  (35.6  meter  1  drms) .  A 
mean  difference  or  offset  between  the  measured  time  differ- 
ence and  the  calculated  time  difference  was  34  nanoseconds 
for    the   9940-X   rate   and   one    nanosecond    for    the   9940-Y    rate. 
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He  also  obtained  a  precision  of  38.0  meters  2  drms  (19.0 
meters  1  drms)  in  the  static  mode.  Nelson  also  states,  that 
the  above  precision  is  only  achievable  if  the  user  has  a 
LORAN-C  receiver  which  has  the  performance  capabilities  of 
those   used      in   the   experiment.  rh9    LORAN-C      receiver   must 

have  "comparable  signal  averaging  time,  extra  notch  filters, 
and    attenua-ion    of    the  signal"   [General    Electric    Co.,    1979]. 

D.       OBJECTIVES 

The  National  Ocean  Survey  requires  that  hyperbolic 
control  systems  used  for  hydrographic  surveying  exhibit  a  1 
drms  of  less  than  0-5  millimeter  at  the  scale  of  the  survey 
[Umbach,  1976].  Although  this  requirement  is  generally  for 
2  mHz  phase  comparision  systems,  it  can  be  inferred  that  it 
also  applies  to  other  hyperblic  systems  such  as  LORAN-C. 
The  scale  routinely  used" for  coastal  surveys  is  between 
1:40,000  and  1:80,000  [Umbach,  1976]  yielding  an  allowable 
error  of  20  to  40  meters  not  including  systematic  errors. 
Schnebele  [1979]  concluded  that  hyperbolic  LORAN-C,  after 
applying  a  single  ASF  Corrector,  is  unsuitable  for  basic 
hydrographic   surveying. 

Whether  or  not  applying  multiple  Additional  Secondary 
Factors  (AS?)  to  LORAN-C  lines  of  position  will  reduce  the 
drms  sufficiently  to  meet  the  accuracy  standards  set  by  the 
National  Ocean  Survey  Hydrographic  llanual  will  be  ascer- 
tained in  this  study.  The  term  multiple  ASF  Correctors 
refers  to  the  applicatioQ  of  more  than  one  corrector  to 
LORAN-C  lines  of  position  over  a  given  area.  The  variable 
ASF  Correctors  result  from  varying  delays  of  the  electromag- 
netic   wave    as   it    propagates    over    different    land   segments. 

Three  methods  of  applying  multiple  ASF  Correctors  were 
tested.      The   first    method   was   the   application    of    a 
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Semi-Empirical  Time  Difference  Grid  Calibration  Model  devel- 
oped by  The  Analytic  Science  Corporation  [1979].  The  spon- 
soring agency  was  the  Onited  States  Coast  Guard.  The 
Semi- Empirical  Model  applies  Secondary  Factors  and 
Additional  Secondary  Factors  for  each  geodetic  position 
based  on  the  distance  over  land,  the  distance  over  water, 
and  the  total  distance  using  mean  se^  water  and  land  conduc- 
tivities. 

The  second  method  which  was  investigated  applies  ASF 
Correctors  found  in  the  DMAHTC  LORAN-C  Correction  Tables 
[DMAHTC,  1981]  to  LORAN-C  lines  of  positions.  These  correc- 
tors were  derived  from  the  ground  conductivities  which  have 
been  determined  in  the  field  by  a  Coast  Guard  calibration 
team    [U.S.    Naval   Oceanograp  hie  Office,    1982]. 

Finally,  a  third  method  was  pursued.  ASF  Correctors, 
which  were  determined  by  field  observations,  were  applied  to 
the    LORAN-C      lines    of    position.  These    were     derer mined    by 

computing  the  difference  between  the  observed  LORAN-C  rates 
and  the  expected  time  difference  which  was  calculated  using 
four  lines  of  position  from  a  very  accurate  microwave  posi- 
tioning system.  These  ASF  Correctors  were  determined  at 
five    minute   latitude    and   longitude   intervals. 
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II.     MIHRJ  Ql  IHI  EIQBLEM 

A.       THE    PRINCIPLES    OF    LORAN-C 

To  understand  the  problems  associated  with  LORAN-C  when 
used  during  hydrcgraphic  operations,  one  must  first  under- 
stand its  principles  of  operation.  LORAN-C  is  a  low 
frequency,  pulsed  signal,  hyperbolic,  radio  navigation 
system,  employing  time  difference  measurements  of  signals 
received  by  the  navigator  from  an  laast  three  ground  tran- 
smitting stations  [Speight,  1982].  The  stations  ar9 
comprised  of  a  master  transmitting  station,  two  or  more 
secondary  transmitting  stations  which  are  strategically 
spaced  several  hundred  miles  apart  and,  if  necessary,  a 
System   Area   Moniter    (SAM)    Station   [U.S.    Coast    Guard,     197U]. 

System  Area  Monitor  (SAM)  stations  associated  rfith  each 
LORAN-C  chain  apply  differential-type  corrections  to  the 
rates  in  real-time.  SAM  stations  continuously  monitor  the 
signals  from  all  transmitters  in  the  chain.  If  the  oberved 
time  difference  deviates  by-  more  than  0.05  jusec  from  the 
expected  value,  then  the  appropriate  secondary  adjusts  its 
emission  delay  time  in  order  to  remove  the  error  [Schnebele, 
1 97  9  ]  . 

The  master  and  at  least  two  secondary  stations  are 
located  such  that  the  signals  from  the  transmitting  stations 
can  be  received  throughout  the  desired  coverage  area.  The 
master  station  is  designated  by  the  letter  "M"  and  the 
secondary  stations  or  slave  stations  are  aesignatisd  W,  X, 
Y,    or    Z    [.U.S.    Coast    Guard,     1974]. 

All  stations  transmit  on  the  common  frequency  of  100 
kHz.  Interference  between  transmitters  is  avoided  through 
the    use    of    time    separation    [Poppe,    1982].         After    the    master 
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station  transmits  a  pulse,  each  secondary  station  delays  its 
own  transmission  for  a  fixed  time,  called  the  secondary 
coding  delay.  This  ceding  delay  is  synchronized  through  the 
use  of  cesium  frequency  standards  at  each  station.  The  high 
stability  and  accuracy  of  these  standards  permit  each 
station  to  derive  its  own  time  of  transmission  without 
reference  to  another  station  [DMA,  1977].  Secondary  coding 
delays  are  predetermined  by  system  propagation  times  and 
equipment   characteristics   [  Laurila,    1976]. 

The  pulse  from  the  master  transmitter  is  distinguished 
from  those  of  the  secondaries  through  phase  coding  of  the 
pulses.  Phase  coding  refers  to  the  inversion  of  the  nega- 
tive and  positive  peaks  of  the  sine  wave  comprising  the  100 
kHz  carrier  portion  of  the  pulse.  The  purpose  of  the  phase 
coding   is   twofold: 

"Firs-,  it  permits  automatic  discrimination  between  the 
master  and  the  various  secondary  stations,  thereby 
permitting  all  stations  to  be  identified  by  their  rela- 
tive  timing    with    respect    to   the   master" 

"Second,  the  phase  coding  nrcvides  protection  against 
excessively  long  skywave  delays  which  would  cause  the 
late  arrival  of  the  proceeding  pulse  to  coincide  with 
the  leading  edge  or  groun dwave  portion  of  a  oulse  being 
tracked"    [Poppe,    1982]. 

The  signals  are  received  by  a  mobile  receiver  where  the 
differences  in  time  of  arrival  of  the  master  signal  and 
various  secondary  signals  are  measurad  and  displayed  en  the 
indicator  portion  of  the  LORAN-C  set.  The  accuracy  of  this 
time  difference  is  increased  by  phase  comparision  "of  the 
synchronized  100  kHz  carrier  within  the  master  and  secondary 
pulses"  [Laurila,  1976].  This  measured  time  difference  (TD 
-  in  microseconds)  represents  a  hyperbolic  Line  of  Position 
(LOP)    [U.S.    Coast    Guard,    197a].       The    intersection    of    two   or 
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Figure   2.1 


Hyperbolic   Fix     (From   Coast    Guard  LORAN-C   User 
Handbook,    197U) 


more      LORAN-C   LOP*s      defines   the      posi*:ion      of   the      observer 
(Figure    2.1).         When    plotted   en    a   chart,    the    intersection    of 
the      resultant      hyperbolic      lines        defines      a      geographical 
position  [Speight,    1982]. 
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B.       PHASE    LAG 

In  a  vacuum,  the  velocity  of  radiated  energy  from  an  antenna 
for  LORAN-C  is  29979  2,458  Icm/sec.  Since  radiated  energy 
cannot  be  shown  pic^orially,  the  phase  of  "che  transmitted 
radiations  is  used.  The  lines  of  constant  phase  of  the 
transmitted  radiation  are  shown  in  Figure  2.2  by  the  curved 
lines  labeled  aa' ,  bb' ,  and  other  similar  designa-ions. 
They  define  the  wave  front  as  it  proceeds  outward  from  zh=i 
antenna  in  all  directions.  The  distance  between  each  line 
of  constant  phase  is  one  wavelength  (x)  [Admiralty  Manual, 
1 96  5  ] . 

\  (meters)    =    29  9792.4  58   km/sec    t    frequency    in   kHz 


"1 


Antenna 


A 


< X ^\< X > 


< X > 


I 


J 


Figure   2.2        Transmitted   Radiation 

The    velocity      of    the   radiated      energy    in   air      depends    on 
:empera-ure,      pressure,        humidity,       and      the   narure      of   the 
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surface  over  which  the  transmissions  pass.  The  retardation 
of  a  -cransmitted  wave  is  known  as  phase  lag.  When  low 
frequencies  are  employed,  such  as  LOHAN-C  at  100,  kHz,  the 
effects  of  change  of  temperature,  pressure,  and  humidity  are 
swamped  by  the  effects  caused  by  changes  in  the  na-ure  of 
the  surface  over  which  the  transmissions  are  traveling 
[Admiralty  Manual,  1965].  The  Genaral  Electric  Company, 
TEMPO      division,  conducted      a      LORAN-C        Signal      Analysis 

Experiment  under  the  direction  of  the  U.S.  Coasr  Guard. 
This  experiment  was  conducted  along  the  U.S.  Wesx  Coast. 
The  General  Electric  Company  recorded  an  overall  change  of 
108  nsec  and  116  nsec  time  of  arrival  from  -tha  master  and 
X-secondary  stations  respectively  after  a  storm  [Samaddar, 
1980].  If  ASF  Correctors  are  as  large  as  two  microseconds 
on  the  West  Coast  [DMAHTC,  1981]  then  the  ASF  Corrections 
are  20  times  larger  in  magnitude  than  weather  effects  for 
the    9940   West   Coast   LORAN-C    chain. 

The  change  in  transmission  ratas  or  phase  lag  are  a 
result  of  the  amount  of  energy  transferred  from  -^he  tran- 
smitted radiation.  This  energy  transfer  depends  en  the 
absorption  qualities  (inversely  related  to  conductivity)  of 
the  surface  over  which  they  are  passing  and  their  wavelength 
(or  frequency) .  The  lower  the  conductivity  and  the  longer 
the  wavelength  (or  lower  the  frequency)  the  greater  the 
transfer  of  energy,  and  vice  versa.  Seawater  has  a  rela- 
tively high  conductivity  (5.0  mhos/meter).  Land  has  a  much 
lower  conductivity,  which  varies  from  marsh  (fairly  high)  to 
dry  sand  and  rock  (very  low)  [Admiralty  Manual,  1965].  Two 
excellant  papers  that  discuss  the  electrical  properties  of 
soil  are  those  of  Smith-Rose  [1934]  and  Pressey,  Ashwell, 
and  Fowler  [1956].  Smith-Rose  [1934]  found  that  the  conduc- 
tivities for  soil  ranged  from  0.18  mhos/meter  for  a  grey 
clay    with   salt   to    .00001    mhos/meter    for    granite. 
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Figure  2.3        Phasa    Lag 

Phase   lag   is      illustrated   in   Figure    2.3,  The   lines   of 

constant  phase  ,  aa',  bb ' ,  and  cc'  become  distorted  as  they 
progress  along  the  sea  surface.  The  dotted  lines  represent 
the  position  of  the  lines  of  constant  phase  in  the  absence 
of  the  sea  surface.  The  wavelength  (X)  in  meters,  measured 
at  heights  of  several  wavelengths  above  the  sea,  remains 
about  the  same  as  the  direct  wave  through  the  earth's  atmo- 
sphere at  299691  km/sec  i-  freguency  in  kHz.  Near  sea  level 
the  absorption  of  energy  retards  -he  progress  of  -he  wave- 
front,  and  makes  the  sea  level  wavelength  (  x' )  less  than  A. 
As  the  lines  of  constant  phase  progress  away  from,  the 
antenna  the  phase  lag  increases  with  distance.  This  is 
known   as   the    Secondary   Factor    (SF)    [Admiralty   Manual,    1965]. 

The  most  intriguing  feature  of  phase  lag  occurs  at  the 
coastal  interface  where  there  is  an  extreme  change  in  the 
conductivites  between  land  and  sea.  Visualizing  the  wave 
front  in  three  dimensions,  the  lower  part  of  the  wave, 
slowed  by  the  drag  of  the  ground,  lags  further  and  further 
behind   the    upper   part    as   the    wave   crosses    the   land.         At    the 
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coastline  it  suddenly  encounters  the  much  lower  impedance  of 
the  sea,  and  in  a  very  short  distance  the  bottom  of  the  wave 
tries  to  catch  up  with  the  top,  as  though  the  whole  wave 
front  were  an  elastic  balloon.  This  is  known  as  "phase 
recovery"  [Eaton,  1979].  Phase  recovery  was  verified  during 
tests  on  Decca  transmissions  across  the  south  coast  of 
England    by    Pressey,    Ashwell  ,    and   Fowler    [1956]- 

The  determination  of  the  Secondary  Factor  for  seawater 
is  fairly  direct  since  the  conductivity  of  seawater  (5.0 
mhos/meter)  is  fairly  constant.  But  for  land  the  conduc- 
tivity can  vary  depending  on  the  type  of  soil  and  its  water 
content    [Smith-Rose,     1934]. 

Phase  lag  for  radiated  energy  over  land  can  be  deter- 
mined   two   ways: 

1)  Assign  an  average  land  conductivity  to  the  ASF 
Model,  For  example,  the  average  conductivity  for 
the  soil  on  the  west  coast  is  0.003  mhos/meter. 
The  average  land  conductivity  will  determine  the 
average  phase  retardation  of  the  path  [The 
Analytic   Science   Corporation    ,1979]. 

2)  Determine  every  conductivity  for  each  portion 
of  a  line  segment  from  the  transmitter  to  the 
receiver.  The  total  of  these  conductivity 
segments  constituting  a  land-water  profile  will 
determine  the  total  phase  retardation  of  the  path 
[Speignt,    1982]. 

C.       TD    MODEL 

Positional  fix  accuracy  using  LORAN-C  is  primarily 
dependent  on  a  chart  makers  ability  to  accurately  compute 
the  expected  difference  in  time-of-acri val  (TOA)  of  received 
groundwave    signals      from    the     transmitting    s-^ations.  Time 
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differences    (TD) ,      are     the    differences   between   the      TO As   of 

the   secondary   and    master   transmitters. 

TD         =  TOA  -  TOA  (2.1) 

i  1  m 

i   =   Secondary    Station 
m   =   Master   Station 

TOA    computations   are   dependent,   upon   an    accurate   knowledge   of 
the    signal    phase   delay. 

The      phase   delay     of   a      groundwave      signal   is      generally 


expressed   as: 


T  +  SF  (2.2) 

nR 


C 


SF 


where  n  is  the  surface  refractive  index,  C  is  the  speed  of 
light  in  a  vacuum,  R  is  the  range  between  the  raceiver  posi- 
tion and  the  trajismi-t;  ting  station.  The  primary  phase  delay, 
T,  is  the  computed  travel  time  of  the  LOEAN-C  pulse  over  a 
distance  equal  to  the  tran smit ter-tD-receiver  greax  circle 
path  length,  taking  into  account  the  veiccitv  of  electromag- 
netic waves  and  the  index  of  refraction  of  the  atmosphere. 
The  secondary  factor  (SF)  is  a  correction  to  the  primary 
phase  delay  and  accounts  for  the  phase  lag.  The  dominant 
term    in    (2.2)         is    the   primary    phase    dslay    (T)  .  The    SF    is 

usually  an  order  of  magnitude  smaller  [The  Analytic  Science 
Corporation,    1979]. 

Thus,   time-of-arr ivals    can   be    expressed   as: 

TOA  =  T  +  SF         +         CD  (2.3) 

i  i  i  i 

TOA  =  T  +  SF  {2,'i) 

m  mm 

where  CD  is  the  true  emission  delay  or  coding  delay  for  the 
LORAN-C  chain  [The  Analytic  Science  Corporation,  1979]-  The 
coding   delay      is    equal  to    a      time   delay    plus   a      computed   one 
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way  baseline  time  (Be)  which  includes  the  secondary  phase 
correction  for  an  all  seawater  path.  The  oneway  baseline 
time  (Be)  is  equal  to  the  distance  between  the  masx.er  and 
secondary  transmitters  in  meters  divided  by  the  propagation 
velocity  of  LORAN-C  through  the  earth's  atmosphere  (299.691 
meters  per  microsecond  [Navigation  Department  DMA,  1982]). 
See    Table      II    for      Coding   Delay      values    for      the    9940      chain 

TABLE    II 
Coding    Delay    9940    LORAN-C    Chain 


Pair  9940-W: 
9940-X: 
9940-Y: 

CD 
CD 
CD 

+  Be  =  1  1000 
+  Be  =  2  7000 
+  Be  =  4  0000 

■••  2796.90  =  13796.90 
+  1094.49  =  28094.49 
+  1967.27  =  41967,27 

|isec 
/isec 
;us6c 

1 

1 

1 
J 

[Riordan,  1979].  Combining  equations  2.1,  2.3,  and  2.4,  the 
true  TD  is  given  by  equation  2.5  [The  Analytic  Science 
Cor  po  rar  io  n ,    1979]. 


TD      =        (T       -      T    )      +     (SF      -       SF    )     +    CD  (2.5) 

i  i  m  i  m  i 


D.       SEMI-SaPIRICAL    TD    GRID    CALIBRATION    MODEL 

The      Semi-Empirical      TD      Model        was      developed      by      The 
Analytic  Science        Corporation  [1979]        in  Reading, 

Massachusetts  for  the  »est  Coast  9940  LORAN-C  chain. 
Similar  "time  difference  (TD)  grid  calibration  techniques 
have  been  successfully  employed  to  develop  an  accurate 
(approximately  100  nsec  drms)  calibra-ed  grid  for  St.  Marys 
River  LORAN-C  chain",  [The  Analytic  Science  Corporation, 
197  9]. 
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^ •      l^chnical    Approach 

The  Semi-Empirical  Model  is  based  on  Millington* s 
empirical  approach  for  computing  the  secondary  factor  over  a 
mixed  (multiple-homogeneous  segment)  path  which  combines 
land  and  sea  phase  delays.  The  generalized  semi-empirical 
polynomial  functional  form  for  the  SF  of  the  LORAN-C  station 
is   given   by: 

K 

2  L 

SF       =    SF(T    ,6)=      I      AT  +S(C      5inl6      ♦    D      ccsl  B  )     (2.6) 

j                    j       j         k=-K^    k    j  1=1       jl               j            jl               j 

where 

j  =  secondary  (W,X,  or  Y)  cr  master  (M)  station, 

nS  , 

T    =   =    jth   station-to   user    primary    phase    delay, 

j        c 

R    =    jth   station-to-user   great-circle   path    length, 
j 

8   =    user    path   bearing   angl=^    at    the    jth   station, 
j 

K    ,    K      and   L    are   oositive   integers, 
12 

C         and   D        are   the   station-dependent   coefficients 
jl  jl 

of   harmonic   terms    in    the   model, 
A      is   the   range-da pandent    coefficient   of    -he    model 

which   may    in   general   be   station-dependent. 

Data  from  27  coastal  sites  distributed  along  the 
West  Coast  and  122  land-sea  sitas  distributed  in  the 
Southern  California  CCZ  (between  Point  Arguallo  and  San 
Diego  -  see  Figure  2.4)  were  used  in  a  Kalman  estimation 
algorithm  to  compute  the  uncertain  coefficients  of  the  land 
and    sea   models      of    the  TD   grid      calibration    algorithm.  (An 

explanation   of   Kalman    filtering    for      the    layman    is   presented 
by   Roger    M.      du    Plessis      [1967].)         The    calibrated   algorithm 
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Figaro    2.4        Loca-iOx'    of    West   Coast    LORAN-C   Statio: 
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was  used  to  compute  TDs  at  each  data  site  and  the  TD  resi- 
duals (difference  between  measured  aad  calibrated  TDs)  were 
examined.  Adjustments  were  then  made  to  the  TD  model  struc- 
ture in  an  attempt  to  further  reduce  the  residuals.  This 
process  of  adjusting  the  aodel  structure  is  repeated  until 
the  residuals  agree  with  the  expected  theoretical  covariance 
associated  with  the  TD  model.  The  model  which  exhibited  the 
"best"  performance  was  selected  as  the  West  Coast  TD  grid 
model. 

2.      Generalized   Ranqe/B earing   Model 

The  Generalized  Range/Bearing  (GRB)  Model  was 
selected  as  the  "best"  semi-empirical  calibration  model  for 
the    West   Coast   chain.      The    semi-empirical    function   is: 

L 
S?      =    AO    +      A1T      +   E     (C      sinlS      +    D      cosl 6  )     usee  (2.7) 

J  j      1=1       jl  j  jl  j      "^ 

where  AO ,  A1 ,  C  and  D  are  the  model  coefficients,  is  the 
path  bearing  angle  measurad  positive  clockwise  from  north  at 
the  jth  (  W,X,Y  or  M)  station  and  I  is  nhe  path  range  to  -he 
jxh    station.  The   GRB    model      is   rela-cively    complex      and    is 

expected  to  exhibit  superior  performance.  The  extensive 
mcdel  is  based  on  knowing  the  distance  overland  (TL)  ,  -he 
distance  over  water  (TS)  ,  the  total  distance  (T)  ,  and  the 
path   bearing  angle    3    (Figtirre   2.5). 

It  was  noted  that  the  calibrated  mcdel  was  expected 
to  be  accurate  and  applicable  only  over  the  extent  of  ranges 
and  bearing  angles  embodied  in  the  calibration  data.  Hence, 
outside  the  region  covered  by  the  calibration  data  the  model 
may  not  be  as  accurate  as  within  the  data  coverage  region. 
Osing  the  GRB  model,  a  drms  value  of  approximately  0.3  jisec 
was  expected  in  areas  where  land  data  alone  was  used  to 
calibrate      the   model.  Inclusion   of      sea    calibration      data 

produced   a    drms    value    of    9.35   to    0.50   jis<5c. 
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The      time      difference      (TD)  is      expressed      by      -he 

following   equation: 

TD      =     (T      -    T    )     +     (SF.     -    SF    )     +    CD.     +    b.         /isec  (2.8) 

i  i  m  11  11 
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where 

n    R. 

T      = |isec, 

i  c 

n    R 
m 

T      = )isec, 

m  c 

R      =   ith   secondary  s-at ion-to-user   great-circle   path   length, 
i 

c      =   -he  speed    of    an   al ectromagnetic   wave   in    a   vacuum 
=    2.99792458    X    10^    m/sec, 

n      =   surface  refractive    index 
=    1.000338, 

R      =   master   station-to- user    great-circle   path    length, 
m 

CD.=  coding    delay    found    in   Table    II, 

b      =   TD   bias  associated    with   the   ith   secondary 
i 

station    (psec)      (Table   III)  . 


TABLE    III 
TD   Bias    (b)     -  jasec 


r- 

"T 

TDW 

-0. 

BbU 

1 

TDX 

-1. 

,173 

1 

c- 

TDY 

-0. 

353 

1 
-J 

SF,=    0.5     {-SI    +    S2    +    S3    -    34    +    S5    +    S6)  (2.9) 


1 


The   term    SI    is   -he    SF    of   a   land    path   of    length    T 
(;isec)     from    the    jth  station: 

SI      =    0.795/T         +    0.439    +     (0.00245)     T  usee.  (2.10) 

sj  sj 

The    terms   S2,      S3,      and   S4      combine   to    make    up   the 
secondary    factior      for    the      seawater    path      Isngihs. 
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S2  is  the  SF  for  the  total  path.  S3  is  the  SF 
using  the  seawater  coefficients  for  the  portion 
with  seawater,  and  S4  is  the  SF  using  the  seawater 
coefficients  for  the  land  path  distance.  The 
seawater  coefficients  are  found  in  Table  IV. 

TABLE  IV 
Seawater  Coefficients 


a1 

a2 
a3 

b1 
b2 

b3 


128.8 
0.  187 
0.000652 

3.188  I 

-0.594  I 

0.000329     I 


SF    (T) 


SF     (T) 

s 


bl/T    +    b2    +    (b3)    T  ;isec, 

if    10    <    T   <   540  ;isec. 


al/T   +    a2    ^-    (a3)    T 
if   T    >    540  jasec. 


;isec. 


(2.1  la) 


(2. lib) 


Term  S5  and  S6  are  the  SFs  of  land  paths  of  length 

T   and  T   .   SF5  is  the  Secondary  Factor  for  the 

J      Ij 
total  length  using  the  land  coefficients  whereas 

SF6  is  the  Secondary  Factor  for  the  distance  over 

land  using  •^.  he  land  coefficients.   The  land 

coefficients  are  found  in  Table  V. 

S5  =  SF  (T  ,  8  )  , 

L   J   J 


S6  =  SF  (T   ,6  )  , 

L   Lj    y 


w  he  r  e 


SF  (T  ,  3  )  =    \0    ^    (A1)  r   +   I  (C   sinl  3  +  D   cosl  3  )  (2.12) 

L   j   j  j    1=1   jl      j     jl      j 
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TABLE    V 


Land  Coefficients 


AO 

= 

1.428 

D 
x2 

= 

0.9U2 

A1 

= 

0.00158 

C 

= 

0.0 

C 

= 

0.0 

c 

= 

0.588 

Wl 

y2 

C 

= 

-0.711 

D 

= 

0.0 

w2 

yi 

D 

= 

0.  323 

D 

= 

0.0 

wl 

y2 

D 

t= 

0.0 

c 

= 

1.010 

w2 

ml 

C 

= 

0.0 

c 

= 

-0.196 

xl 

in2 

C 

= 

0.0 

D 

= 

-0.893 

x2 

ml 

D 

= 

0,0 

D 

= 

-0.355 

Xl 

m2 

._    J 

S.       DMAHTC    MODEL 

1 .  S6 a   S F    Model 

The   equations    for   the    Sea   SF    Model    is: 

SF    =     (B1/T)     +    32    +     (B3   T)    usec,    if    ID    <    T    <   537  jis^c,     (2.13a) 
SF    =     (A1/T)    +    A2    +     (A3   T)     a  sec,    if   T    >    537  jtisec,  (2.13b) 

where  T  is  -he  primary  phase  delay  (or  range)  in  microse- 
conds (usee);  Ak  and  Bk  (k  =  1,  2,  and  3)  are  the  sea  model 
coefficients  used  by  DMAHTC  in  program  TDGRID  [Funakoshi, 
1982].      The  coefficients   are   found   in    Table   VI. 

2.  Land    SF    Model 

The      solution      used      to    resolve      DMAHTC      Taoles      ASF 
Corrections   is      called  Millington's    Method      [DMAHTC,       1981]. 
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TABLE    VI 


Sea   SF   Model  Coefficients 


A1 
A2 
A3 

B1 
iB2 

|B3 


129.04323 
-0.40758 
0.00064576813 

2.741282 
-0.011402 
0.00032774815 


— I 


This  method  is  based  on  the  premise  that  the  phase  distor- 
tion due  to  a  composite  land-sea  path  is  the  arithmetic 
average  of  the  phase  distortion  foand  in  the  forward  and 
reverse  paths  of  the  propagated  signal  [  DMAHTC,  1981].  For 
example,  in  Figure  2.6  two  azimuths  have  been  drawn  on  the 
map  and  are  labeled  as  210°  and  235o.  Also  placed  on  the 
map  are  the  proper  ground  conductivities  which  have  been 
determined  in  the  field  by  the  Coast  Guard  calibration  team. 
A  great  circle  drawn  on  the  appropriate  chart  or  charts  from 
the  LORAN-C  Station  coordinates  to  the  area  undsr  considera- 
tion spans  various  lengths  of  land  and  seawater.  Each 
length  or  segment  will  have  a  specific  conductivity  and 
distance.  The  total  of  these  conductivity  segments,  consti- 
tuting a  land-water  profile,  will  determine  the  total  phase 
retardation  along  that  path  [U.S.  Naval  Oceanographic 
Office,  1982].  All  azimuths  and  distances  are  computed 
based  on  the  World  Geodetic  Systsm  (WGS)  datum  [DMAHTC, 
198  1].  The  values  of  phase  retardation  for  a  given  ground 
conductivity  are  tabulated  in  the  National  Bureau  of 
Standards    (N3S)     Circular    57  3   [Speight,    1982]. 

The  formula  used  to  derive  the  ASF  Correction  for 
the  time  difference  for  a  master-slave  transmitting  station 
pair    is: 

ASF   Correction   =    (-Slave   Error)    -    (-Master   Error) 
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TRANSMITTER 


(,003  mhos/m) 


SEAWATER 
(5.0   mhos/m) 


^ 


Figure    2.6        Composite   Land-Sea    Path 

The   mean      values    derived      for    one      station    from  the 

forward     and    reverse      solution      of      Millington's    Me-hod  are 

subtracted      from   the      Sea  SF      Model.  The    differences  are 
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presented  in  an  azimuth  array.  This  array  is  a  series  of 
geodetic  azimuths  radiating  from  the  transmitter  with 
corrections  computed  at  incremented  distances  along  each 
azimuth.  After  the  valaes  are  computed  in  the  azimuthal 
array  a  compu-cer  program  rearranges  them  into  a  matrix  form. 
The  matrix  form  is  the  arrangement  of  corrections  into  rows 
and  columns  covering  a  specified  geographical  area  at  a 
constant  spacing.  The  purpose  of  the  matrix  is  to  enable 
the  corrections  from  two  LORAN-C  transmitters  to  be  added 
algebraically,  combined  into  a  single  matrix,  and  arranged 
in  the  desired  tabular  form.  This  tabular  form  is  the  body 
of    the   table   [DMAHTC,     1981]. 

3 •      Ta b le   Description 

Each  table  contains  a  complete  chain.  Figures  2.5, 
2.6,  and  2.7  depict  LORAN-C  ASF  Correctors  for  chains 
9940-W,  X,  and  Y  for  Monterey  Bay,  California.  A  table 
section  is  prepared  for  each  station  pair  (master  station 
and  one  slave  station)  in  a  LORAN-::  chain.  As  a  rule  the 
limits  of  the  table  coverage  are  determined  by  the  range  of 
the  groundwave  transmissions  for  t-hs  LORAN-C  chain.  Each 
page  of  corrections  in  the  table  covers  an  area  three 
degrees  in  latitude  by  one  degree  of  longitude,  with  correc- 
tions printed  in  increments  of  five  minutes  of  arc.  Rate 
designation  and  page  numbers  are  printed  at  the  top  of  each 
correction  page.  Those  pages  where  latitude  and  longitude 
limits  contain  both  land  and  sea  are  included  but  correc- 
tions apply  only  for  the  area  covered  by  the  a.S.  Coastal 
Confluence  Zone  (CCZ) .  Large  land  bodies  and  areas  ouside 
the  CCZ  are  represented  by  blank  spaces  on  the  page.  ASF 
correcticn  values  can  be  either  positive  or  negative  (posi- 
tive values  are  shown  without  sign).  Areas  requiring  no 
correction    show    a    zero   value    which   in    some    cases    in    preceded 
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Figure    2.7 


9940-v*    ASF   Correction   Taole 


by  a  negative  sign  indicating  that  the  zero  results  from  the 
rounding  off  of  a  value  slightly  less  than  zero  (indicates 
the   trend   of   the   correction). 
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Figure   2.8        99aO-X  ASF  Correction   Tacls 

The  table  can  be  entered  directly  by  using  the 
ship's  position  determined  to  the  nearest  five  minutes  of 
arc    in      la-irude   and   longitude      either    by    dead     reckoning    or 
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Figure    2.9        9940-Y   ASF   Correction   Table 

scr^e  other  means.  To  find  the  page  with  the  appropriate 
correction,  the  Page  Indexes  of  the  table  should  be 
utilized.        These    indexes   show   the    limi-s    and   page   number    of 
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all  pages  in  the  table.  To  locate  the  number  of  the  page  on 
which  the  desired  correction  is  to  be  found  ihe  Page  Index 
is    entered      with   the      ship's   position.  In   some      cases   the 

ship's  position  will  fall  on  the  page  limit  in  either  lati- 
tude or  longitude  or  both.  These  positions  are  repeated  on 
both    pages    and   either    page    may   be  used. 

The  ASF  Correction  is  added  algebraically  to  the 
time  difference  for  the  LOR AN-C  pair.  Interpolation  of  this 
data  will  not  necessarily  improve  the  accuracy  due  to  the 
method  used  to  determine  ASF  Correctors  [DMAHTC,  1981]. 
Since  the  correctors  are  computed  in  the  azimuthal  array  and 
are  based  on  the  conductivity  and  distance  over  which  the 
LORAN-C  electromagnetic  wave  travels,  the  ASF  Corrector 
between  the  published  AS?  Correctors  in  the  tables  may  not 
be  the  linear  interpolated  values.  For  example,  the  ASF 
Corrector  for  a  distance  of  500  m  with  an  azimuth  of  180°  is 
equal  to  1.5  jisec.  The  ASF  corrector  for  a  distance  of  500 
m  with  an  azimuth  of  181^  is  equal  to  1.6  psec.  The  inter- 
polated value  between  1.5  and  1.6  is  1.55.  The  true  ASF 
Corrector  for  the  180°  30*  azimuth  is  1.U  since  the  land 
distance  for  the  same  azimuth  is  less  than  the  land  distance 
for  the  180O  and  181 o  azimuth.  The  LORAN-C  signal  passed 
over    a   harbor   [Dansford,    1982]- 

F.  ATTEMPTED  DETERMINATION  OF  ASF  CORRECTORS  BY  HYDRO  FIELD 
PARTIES 

One  of  the  major  problems  encoantered  by  hydrographic 
survey  operating  units  when  using  LORAN-C  for  position 
control  is  the  determination  of  the  AS?  Correctors  for  the 
survey  area.  The  four  surveys  mentioned  in  the  Introduction 
all  made  attempts  to  determine  the  correctors  by  comparing 
the    LORAN-C  rates    to    a  second  source. 
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Bathymetric  Surveys  H-9822  [NOAA  H-9822,  1979]  and 
H-9823  [NOAA  H-9823,  1979]  Gulf  of  Alaska,  compared  the 
rates  from  an  Internav  LC-20U  LORAN  receiver  to  computed 
rates  from  a  position  obtained  from  a  JMR-1  Satellite 
Navigation  Receiver  when  available.  Shore  ties  using  radar 
ranges,  visual  bearings,  and  sextant  angles  in  comparison  to 
LORAN-C  rates  were  also  made  prior  to  and  after  each  survey. 
The  calibrations  of  LORAN-C  rates  were  based  on  the  satel- 
lite positions  only  since  the  positional  computation  of 
LORAN-C  and  JMR  Doppler  Satellite  were  made  on  the  WGS  1972 
datum  whereas  the  land  ties  were  based  on  the  NAD  1927 
datum. 

Bathymetric  survey  SU- 40-7-82  which  extended  along  the 
Washington,  Oregon,  and  California  Coasts  used  LORAN-C  as 
navigational     control.  LORAN-C      time        differences      were 

compared  with  SATNAV  positions.  The  report  did  not  indicate 
whether  any  correctors  were  applied  [NOAA  Ship  Surveyor, 
198  2]- 

Finally,  Hydrographic  Survey  S-K902-Wh-82, 

Reconnaissance  Survey  of  Safety  Fairways,  Gulf  of  Mexico 
used    LORAN-C     as   a    pcsizioning      control    system.  The    posi- 

tioning unit  was  an  LC-20'i  receiver.  LORAN-C  rates  were 
input  via  the  HIDROPLOT  Controller,  a  special  purpose 
input-outpu"*-  interface  which  is  the  nucleus  of  the  computer 
system  hardware  [dmbach,  19  76].  Positions  were  computed  and 
plotted  by  Program  RK121,  LORAN-C  Real-Time  HYDROPLOT 
[Backus,    1980]. 

ASF  Correctors  for  LORAN-C  were  achieved  by  visual  cali- 
bration using  three  poirn:  sextant  fixes  using  charted  oil 
rigs   as     control   in      the   vicinity      of    the      survey    area.  A 

three-point  sextant  fix  is  a  convenient  and  accurate  method 
for  determining  the  position  of  a  hydrographic  survey 
vessel.         Sextants      are   used   to      measure    two      angles    between 
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three  objects  of  known  geographic  positon.  The  center 
object  is  common  to  both  angles.  The  position  of  the  obser- 
vers taking  rhe  angles  is  fixed  by  the  intersection  of  three 
circular   lines  of    position    [Umbach,    1976]. 

These  sextant  angles  were  recorded  and  later  tranf erred 
onto  their  respective  charts  using  a  plastic  rhree-arm 
protractor.  A  plastic  three-arm  protractor  is  transparent 
and  made  up  of  one  fixed  arm  and  two  movable  arms  which 
contain  an  etched  line  that  is  radial  with  the  center  of  the 
protractor  [Ombach,  1976].  Sextant  angles  observed  in  pairs 
for  a  resection  fix  with  a  common  center  mark  may  be  plotted 
directly  by  this  instrument.  When  the  three  arms  are  placed 
at  the  angles  observed  and  fitted  sd  as  to  pass  through  the 
plotted  positions  of  the  observed  stations  on  the  field 
sheet,  the  hole  at  the  center  of  the  three-arm  protractor  is 
the    fixed    positon   of    the   vessel  [Ingham,    1975]. 

Partial  correctors  for  each  area  surveyed  were  defined 
by  comparing  the  observed  rates  and  the  determined  rates 
plotted  on  the  nautical  chart.  The  partial  correctors  were 
applied  via  the  HYDROPLOT  Controller.  However,  even  after 
applying  these  correctors,  the  plotted  position  still  disa- 
greed with  the  ship's  determined  postion  with  respect  to  the 
oil  rigs.  Ship's  personnel  attributed  the  discrepancies  to 
one    or   more   of   the    following: 

1)  Accuracy   of   the   charted   rigs, 

2)  Weather    effect    on   LORAN-C, 

3)  Time    of    day, 

4)  Propagation    of    signal   over    land   path, 

5)  Three-arm   protractor   accuracy,    and 

6)  Error    in    the   conversion    by   the   software 
of    the   LORAN-C    rates   to    latitude   and 
longitude   [NO\A,    1982]. 
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There  is  an  apparent  need  for  a  LORAN-C  calibration 
routine  aboard  NOAA  ships  which  provides  the  ASF  Correctors 
for  program  RK121,  LORAN-C  Real-Time  HYD20PL0T.  The  routine 
should  use  the  same  geodetic  distance  computation  found  in 
RK121  and  use  the  same  datum  as  that  of  the  nautical  chart 
of   the   survey  area. 

The  above  mentioned  discrepancies  illustrate  the  defi- 
ciencies in  applying  a  single  ASF  Corrector  -o  LORAN-C  data. 
The  accuracies  for  hydrography  cannot  be  met  using  single 
correctors  because  the  errors  are  non-linear  and  systematic. 
They  cannot  be  distributed  like  residuals  in  a  traverse. 
Schnebele  [1979]  has  already  proven  that  single  ASF  area 
correctors  to  LORAN-C  positions  do  not  meet  the  accuracy 
s-andards   of   the   NOS    Hydrographic  Manual. 

Based  on  visual  inspection  of  the  DMAHTC  LORAN-C 
Correction  Tables,  ASF  Correctors  should  be  updated  every 
five  minutes  of  Latitude  or  Longitude  change.  In  Monterey 
Bay,  California,  there  is  approximately  0.1  to  0.2  usee 
difference  for  every  five  minures  of  change,  a  poten-ial 
error   of   55  to    110    meters. 
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III.       EXPERIMENIAL    PROCEDURE 

A.       FIELD    PROCEDURES 

In  order  to  compare  the  use  of  differential  LOHAN-C  with 
ASF  multiple  correctors,  typical  survey  operations  were 
planned  for  the  southern  portion  of  Monterey  3ay, 
California.  This  survey  was  conducted  in  con juncxion  with  a 
comparative  evaluation  of  multiple  lines  of  positon  for 
selective  positioning  methods  [Anderson,  1982].  Four  micro- 
wave  ranging   systems   were   set      on   Known   geographic    positions 

TABLE   VII 
Geographic   Names    and    Positions 


Microwave    System 

Stations 

Seaside   4     (1964) 


Use    Hon   Ecc. 
Geoceiver    Sec. 
Mussel   Sec. 
Park    (1931) 
Mulligan   RM1 
Range   7    (197  2) 
Mussel     (1932) 


Geographic   Position 
(NAD   19  27) 


35  0  36 
1210  51 

360  36 

1210  52 

36  0  36 
121  0  53 

36  0  36 

1210  54 

360  53 

121  0  49 

360  44 

121  0  47 

360  39 

1210  49 

36  0  37 

1210  54 


23.44596" 
38.83281" 

04,73031" 

35.98040" 

32.49281" 
25.21162" 

18.25484" 
11.49661" 

13.30600" 

46.743  00" 

5  6.49531" 
52.31090" 

02.47787" 

08.58202" 

18.15100" 
11.49661" 


^§l£§   !^§^d. 
June   3-5,1982 

June    3-5,1982 

June    3-5,1982 

June    3-5,1982 

June    6-7,1982 

Jane    6-7,1982 

June    6-7,  1982 

June   6-7, 1982 
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listed  in  Table  VII.  A  series  of  tracklines  were  run  in  two 
separate  areas  as  shown  on  ?igure  3-1.  To  ensure  that  the 
microwave  postioning  system  was  working  properly,  the  equip- 
ment was  calibrated  over  known  baselines  of  1497.47  meters 
and  7877.3  1  meters  at  the  beginning  and  end  of  the  project. 
Trackline  observations  were  only  made  during  the  daytime  in 
fair  weather  conditions  so  as  to  eliminate  sky  waves  and 
weather  changes  that  influence  LORAN-C  signal  propagation 
characteristics  [Samaddar,  1982].  The  vessel  used  was  the 
126  foot  R/V  Acania  which  is  operated  by  the  Naval 
Postgraduate   School. 

The  positioning  equipment  consisted  of  a  Micrologic 
ML- 1000  LORAN-C  receiver  (0.0  1  usee  resolution)  and  a 
Trisponder  Microwave  System  provided  by  Racal-DECCA  Survey, 
Inc.  The  Trisponder  Microwave  System  consisted  of  four  DNT1 
Model  217C  transponders,  four  DNT1  Model  21017  HP  sector 
antennas  with  87o  by  5°  beam  widths,  cne  DNT1  DDMU  (Digital 
Distance  Measuring  Unit),  two  Omni  DVTI  Model  2  1019  HP 
antennas,  a  Houston  Instruments  Model  DP3-M2D/RC3  plotter, 
and   a      Texas   Instruments      743   terminal       (Table   7III) .  The 

manufacturer's  published  accuracy  for  the  positioning  equip- 
ment is  ±  1  m  for  a  single  range  [  Ra::al-DECCA  Survey,  1981]. 
Anderson  [1982]  discusses  the  accuracy  of  four  lines  of 
position.  The  four  Decca  Trisponder  distances  were  recorded 
via  a  Texas  Instruments  743  data  terminal  while  the  LORAN-C 
rates  were  manually  logged.  The  data  was  acquired  at  one 
minute  intervals  while  the  ship  maintained  constant  course 
and  speed.  The  recorded  LORAN-C  rates  were  9940-W,  X,  and  Y 
of   the   West   Coast   chain. 

To  test  the  potential  for  calibrating  the  LORAN-C  System 
using  the  Semi-Smpirical  Model,  the  correction  tables,  and 
multiple  observed  field  correctors,  the  positions  derived 
from    the   microwave   system    n easurement s      were   used    to    compute 
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Figure    3.1        Location   of  Test   irsas    in   Mcntsrey   Bay 
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TABLE    VIII 


Microwave   Positioning   Equipment 


Eauipment 

S/N 

DNT1    Model   217C   Transponders    (Code 

(Code 
(Code 
(Code 

72R] 
74R 
76R 
78R 

1 
1 

3323 
3320 
3321 
3322 

DNT1    Model   2107    HP    Sector    Antenna 
(870    by    50    Beam    Widths) 

185 
186 
187 
191 

DNT1    DDMV    Model    540 

426 

Omni,    DNT1    Antenna   Model   2  1019    HP 

194 
200 

Houston   Instruments   Plotter   Model 

jP3- 

12D/RC3 

10722-10 

Texas   Instruments    7U 3  Terminal 

34418 
—          _    _          _j 

expected      LORAN-C      time      differences    at      each      point.  The 

difference  or  offset  between  these  expected  time  differences 
and  the  observed  values  were  computed  for  the  three  methods. 
The  mean  offset,  standard  error,  and  drms  values  were  also 
computed  and   compared. 

3.       MICROWAVE    SYSTEM    POSITIONING 

The  geographic  position  of  the  ship  based  on  four  lines 
of  position  was  determined  using  a  computer  program  callei 
GP3YLQ  (Geographic  Position  by  Least  Squares)  written  by  the 
author  (See  Appendix  A)  .  GPBYLQ  contains  subroutine  LSQR 
(Least  Squares),  which  is  a  least  squares  adjustment  written 
by  Paul  R.  Wolf,  Ph.D.  [1974]  and  revised  by  LCDR  D.  Leath 
[1981].  Geographic  postions  were  converted  to  X,Y  (meters) 
which  in  turn  were  converted  to  geographic  position  via 
subroutines   GPTOXY    and  XYTOGP,    respectively   [Wallace,    1974], 
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Subroutine  GPTOXY  and  XYTOGP  are  based  on  the  Modified 
Transverse  Mercator  Grid  srfhich  was  centered  in  the  survey 
area. 

The  Modified  Transverse  Mercator  (MTM)  projection  is 
used  by  the  National  Ocsan  Survey  and  is  similar  to  the 
projection  used  in  the  Universal  Transverse  Msrcaror  (UTM) 
system.  The  main  difference  is  that  in  the  MTM  a  Central 
Meridian  is  picked  that  is  near  the  survey  area  instead  of 
being  fixed  at  a  particular  meridian  [Wallace,  1971]. 
Central  Meridian  (CMER)  ,  False  Easting  (FEST)  ,  and 
Controlling  Latitude  (CLAT)  are  the  three  parameters  which 
define  the  MTM  projection.  CMER  is  the  mean  longitude 
computed  using  the  maximum  and  minimum  longitudes  of  the 
survey  limits,  FEST  is  the  X-Coordinate  that  is  assigned  to 
the  Central  Meridian,  and  CLAT  is  the  distance  in  meters 
from  the  equator  to  seme  reference  latitude  [Wallace,  1971], 
The  Central  Meridian,  False  Easting,  and  Controlling 
Latitude  used  for  Monterey  Bay,  California  referenced  to  NOS 
Chart    18685   are: 

CMER    =    1210    56'     00  .0", 

FEST    =    20000  .0, 

CLAT    =    4050000.0, 
To   be   consistent      with    the   National    Ocean      Survey   charts 
of   the   area,        all    computations    were    done      relative    to    North 
American      Datum       (NAD)         1927      geographic      positions.  All 

programs  were  executed  on  an  IBM  30  33  computer  located  at 
W.R.  Church  Computer  Center,  Naval  Postgraduate  School, 
Monterey,    California. 

C.   LORAN-C  COMPUTATIONS 

The  differences  or  offset  between  the  observed  and 
computed  LORAN-C  rates  using  the  Semi-Empirical  TD  Model, 
ASF      LORAN-C   Correction      Tables,      or      the    Multiple      Observed 
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TABLE    IX 


LORAN-C    9  940   Chain    Data 


r 

Station 

Geogrj 

aphic   Position 
(NAD   1927) 

Master   -   Fallon,    Nevada 

390 
1180 

33' 
49' 

07.03"N 
52.23"W 

Slave    -   George.    Washington 
9946-W 

470 
1190 

03' 
44" 

48.82"N 
34.78"W 

Slave   -   Middletown, 

California   99 40-X 

380 
1220 

46' 
29' 

'    57.49"N 
40.04"W 

Slave    -    Searchlight,    Nevada 
9940-Y 

350 
1140 

19 

48' 

•     1 8 . 3  2  ••  N 

'     13.95»W 

t -                                

_      J 

Correctors  were  compared  to  the  offsets  between  the  observed 
TD  rates  and  the  computed  rates  for  which  only  the  seawater 
Secondary  Factors  (SF)  were  applied.  The  comparison  of  the 
offsets  between  the  four  methods  illustrates  the  improvement 
in   positional   accuracy  after   applying    ASF   Correctors. 

^  •      Sea  water    Secondary    Factors    (SF) 

Time  differences  using  only  the  seawa'rer  Secondary 
Factors  for  each  of  the  geographic  positions  were  computed 
using  program  LORAN  written  by  the  author  {Appendix  C) . 
Seawater  Secondary  Factors  (SF)  were  computed  using  formula 
2.13      and      the   coefficents      found      in      Table    VI.  All      TOA 

distances  in  meters  were  determined  using  subroutine  INVER1. 
INVER1  is  a  geodetic  inverse  routine  using  T.  Vincenty's 
modified  Rainsford's  method  with  Helmert's  elliptical  terms, 
programed  by  LCDR  L.  Pfeifer,  NOAA  [1975].  Subroutine 
INVER1  is  accurate  to  0.0  001  m  halfway  around  the  world 
[Pfeifer,  1982].  All  distances  were  converted  to  microse- 
conds  using    299.792458   m/|iS€C.         Time    differences     (TD)       were 
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computed      from  equation      2.5.  North      American    Datum      1927 

geographic  positions  were  used  for  all  computations.  See 
Table  IX  for  the  positions  of  LORAN-C  9940  transmitters 
[Riordan,    1979]. 

To  ensure  that  subroutine  INVER1  was  functioning 
properly  the  distances  between  ths  master  and  secondary 
stations  were  compared  to  the  NOS  published  baseline 
distance  [Riordan,  1979].  The  published  distances  and  -he 
results   from   routine    INVER1    are    lisred   below. 

Published   Baseline  Commuted    Baseline 

Dis-ance     (m)  Distance    z   INVERT     (m) 

99U0-W  837,777.0929  837,777.115 

99U0-X  327,886.3720  327,886.316 

99aO-Y  589,298.5712  589,298.589 

The  difference  between  the  publishei  and  computed  baselines 
ranged   from   0.02   to    0.06    m. 

Differences  or  offsets  (x  .  )  were  obtained  by 
subtracting  the  observed  LORAN-C  rates  from  the  computed 
values   from      the   various   methods.  The    mean      difference    or 

offset  (X)  and  standard  deviation  (s)  in  microseconds  for 
each  rate  were  determined  using  equations  3.1  and  3,2 
[Wonnacott,    1935]: 

n 
1 
Mean  offset    (/isec)  :  x   =   -      Z      x  (3.1) 

n  i 

i=1 

where:  x   =  original  observation  in  usee, 
i 

n   =  number  of  observations; 


n 
Standard  1 

error:  s   =  Z  (x   -  x)2.      (3.2) 

(;isec)  n-1      i 
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The    mean  offset      in  microseconds   can   be      converted   to   meters 
using   equation   3.3   [Heinzen,    1977]: 


Mean  offset 
(meters) 

where  w 


1  = 


X      w 

sin  01 


(3.3) 


=    distance    corresponding   to   one    microsecond 

on    the    baseline 
=    149.396229   meters    [Bigelow,    1963]. 

=   one-half    the    angle    between   the 

radius   vectors   from    the    position   to 
the  master  and   secondary   stations. 


Equation   3.4    was   used    to    compute   drms    values    in    meters 
[Heinzen,    1977]: 


drms= 


sin  (      +        ) 
1         2 


\1 


(S^W)2 

2 

sin  a 
1 

(S^W)2 

+ + 

2 

sin  a 
2 

2  pcos ( a    +  a     )      s    w 
1         2          1 

\ 

s    w 
2 

(sin  01  )     (sin  a  ) 

(3.4) 


where: 


=  correlation   coefficient   =    0.33, 


s    and    a   are    as   above   with    the    subscript 

denoting    the   appropria-^e 
secondary   station. 

The  correlation  coefficient  (p)  is  a  resuli  of  the  secondary 
station  having  a  common  line  of  position  with  the  master 
station.  Although  often  ignored,  various  authors  assign 
values  ranging  from  0.33  to  0.40.  Bigelow  [1963]  chooses 
P      =    0.33. 

2  .      Semi-Empirical  TD  Grid 

To  determine  if  the  Semi-Empirical  TD  Grid  would 
reduce  either  the  drms  value  or  the  offset  between  the 
observed  and  expected  TD  rates,  program  LOPLC  (Line  of 
Position  -  LORAN-C)  was  written  by  the  author  (Appendix  B) . 
Program  LOPLC  computes  the  distance  over  land,  distance  over 
water,  the  total  distance,  and  the  azimuth  from  north  for 
the      transmitting      station    using      Subroutine      INVER1.  The 
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land-sea  distances  were  computed  by  selecting  a  point  which 
was  located  along  the  coast  of  Monterey  Bay,  Subroutine 
SHORPT  (Shore  Point)  interpolated  a  geodetic  position  from 
23  geodetic  points  which  outlined  Monterey  Bay,  California, 
The  land/sea  Secondary  Factor  was  determined  in  Subroutine 
SECFAC  which  is  found  in  program  LOPLC  using  eguations  2,8 
through  2.12,  All  corapu taxions  ware  based  on  NAD  1927 
geodetic   datum. 

Tests  were  made  on  program  LOPLC  using  data  found  in 
the  Semi-Empirical  TD  Grid  article  [The  Analytic  Science 
Corporation,  1979]  using  the  WGS  1972  datum.  Station  TASC 
55  located  at  latitude  3ao  34*  18,3"  N  and  longitude  120o 
39*       40. 3"      W,  was    selected      from      The      Analytic      Science 

Corpora-ion   article.  It    was  one      of   the   stations      used   to 

calibrate  the  coefficients  for  the  Semi-Empirical  TD  equa- 
tions discussed  earlier.  The  only  distances  listed  for  TASC 
55  were  the  individual  distances  ov=r  land  and  over  water. 
The  total  distance  between  the  transmitters  and  TASC  55  was 
computed     by     adding      the      land      and      sea      distances.  The 

following  station-tc-site  pazh  segment  lengths  are  listed 
for    TASC   55: 

TASC    55  OI2R2 

Distance  Dist ance  Distance  Distance 

Station                    (km)  (km)  (km)  (km) 

Master                    540,730  35.548  575.278  576.083 

X-Secondary        369,248  126,096  495.344  495.178 

Y-Seccndary        525,659  15,832  541.491  541.308 

Using  the  land  distances,  sea  distance,  and  total  distances 
to  TASC  55,  and  the  computed  azimth  from  subroutine  INVER1, 
program  LOPLC  produced  offsets  of  3.86  jisec  and  1.13  jisf^c 
for  9940-X  and  9940-Y  respectively  when  compared  to  the 
expected   time   differences   at   TASC   55. 


50 


An  attempt  was  made  to  determine  if  this  discrepency 
was      due      to   program      LOPLC.  A      comparison   of      the      total 

distances  in  the  report  between  the  transmitters  and  TASC  55 
were  compared  to  the  total  distances  computed  by  INVEE1. 
The  difference  in  distance  between  Subroutine  INVER1  and  the 
total  distances  from  the  Semi-Empirical  TD  3rid  acticle 
ranged  from  160  -  200  m  corresponding  to  0.5  to  0.65  usee. 
This  results  in  time  difference  errors  of  0.04  usee  for  the 
9940-X  rate  and  a  0.10  ;isec  for  the  9940-Y  rate,  implying 
that  the  offsets,  0.86  psec  and  1.13  fisec,  are  caused  in 
part  by  the  method  in  which  the  azimuth  from  north  is  deter- 
mined. Unfortunately,  azimuth  data  from  TASC  55  was  not 
presented  in  the  article.  The  Ar;alytic  Science  Corporation 
has  been  contacted  on  numerous  occasions  in  an  attempt  to 
ascertain  their  method  of  determining  distance  and  azimuth. 
As   of   this    date  there    has   been   no   response. 

Nevertheless,  ttie  data  from  Monterey  Bay  was 
utilized  in  Program  LOPLC  to  obtain  results  that  could  be 
compared   to      that   obtained    by      the   other    methods.  If   this 

method  is  accurate  enough,  AS?  Correctors  could  be  deter- 
mined via  computer  for  each  individual  position  without 
using  tables  or  field  determined  correctors.  The  mean 
offset  (X)  and  standard  arror  (s)  in  microseconds  between 
the  observed  and  calculated  rates  w=re  computed  using  equa- 
tions 3.1  and  3.2,  respectively.  The  mean  offset  (1)  and 
drms  in  meters  were  computed  using  equation  3.3  and  3.4, 
respectively. 

3.      Calculated   Table    ASF   Correctors 

The  offset  between  the  observed  LORAN-C  rates  and 
the  expected  values  with  applied  ASF  Correctors  from  the 
LORAN-C  Correction  Tables  and  the  seawater  Secondary  Factor 
(equation    2.13)       were    determined      using    program    LORTAB    which 
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was  written  by  the  author  (Appendix  D)  .  The  ASF  Corrector 
for  each  position  was  selected  using  subroutine  TABLE. 
Subroutine  TABLE,  which  is  found  in  program  LORTAB  deter- 
mines an  AS?  Corrector  for  each  data  point  based  on  its 
geodetic   position.  The   ASF  Correctors   used      in    subroutine 

TABLE  (see  Figures  2.5,  2.6,  and  2.7)  are  located  between 
latitude  36°  35'  N  and  36o  55'N  and  longitude  121°  50'  W  and 
1220      00'      W.  The      difference   in      sign      between      the      ASF 

Correctors  in  subroutine  TABLE  and  those  found  in  ths 
LORAN-C  Correction  Table  is  due  to  the  difference  in  their 
application.  ASF  Correctors  from  the  tables  are  applied  to 
observed  rates  while  ASF  Correctors  from  Subroutine  TABLE 
are  applied  to  the  calculated  LORAN-C  rates.  Negative  ASF 
Correctors  from  the  LORAN-C  Correction  Table  were  applied  to 
the  calculated  time  differences  ro  be  consistent  with  the 
application  of  Secondary  Factors  to  the  computed  primary 
phase   delay. 

As  before  equations  3.1  through  3.4  were  used  to 
compute  "he  the  mean  offset  ,  standard  error,  and  drms.  If 
this  application  of  LORAN-C  Correction  Tables  is  accurate 
enough,  it  precludes  the  need  to  determine  ASF  Correctors  in 
the    field. 

'*  •      Observed   ASF   Cor rec  tors 

Observed  ASF  Correctors  were  determined  using 
Program  ASFSEL  (ASF  Selection)  which  was  written  by  the 
author  (Appendix  E)  ,  Schnebele's  prior  data  [1979],  and  the 
June  1982  data.  Program  ASFSEL  (ASF  Selec-ion)  was  writ-en 
by  the  author.  This  program  computes  the  ASF  Correctors  by 
subtracting  the  observed  LORAN-C  rates  from  the  expected 
values.  Only  the  seawater  Secondary  Factors  from  equation 
2.13    have   been   applied  to   the  calculated    time    differences. 
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The  mean  ASF  Correctors  for  the  LORAN-C  rates  were  deter- 
mined at  every  minute  of  latitude  and  longitude  between  36° 
50*    N      and    36°      35'    N      and    122°      04»W    and      122o    49'       W.       See 

j  Appendix  E    for   mean   ASF     Correctors   at    one   minute    intervals. 

The  ASF  Correctors  were  then  selected  and  assigned  to 
subroutine  TABLS  in  Program  L0RTA3  at  five  minute  latitude 
and   longitude      intervals.         All   values      were   entered      to   the 

I  nearest    hundredth   of    a      microsecond.         Equations    3.1    through 

3.4  were  used  as  before  for  computations.  This  determina- 
tion of  ASF  Correctors  in  the  field,  if  accurate  enough,  may 
allow  the  use  of  LORAN-C  as  a  positioning  system  for  hydro- 
graphic  surveying   in    the    future. 


53 


IV.       RESULTS 

A  total  of  620  time  differences  (TD)  and  geographic 
positions  based  on  four  lines  of  position  were  recorded  in 
the  southern  portion  of  Montery  Bay  in  order  to  compare  the 
use  of  differential  LORAN-C  with  ASF  Multiple  Correctors. 
The  data  was  divided  into  four  sets.  The  first  data  set  is 
Schnebele's  [1979]  prior  data  consisting  of  a  tozal  of  130 
data  points  collected  on  two  separate  days,  Juna  12  and  July 
25,  1979  between  latitude  36°  38»  N  and  36°  47'  N  and 
longitude  121o  49'  W  and  12 20  02'  W.  The  recorded  LORAN-C 
time   difference   rates    were    9940-Y   and   9940-W. 

The  second  data  set  consists  of  193  time  differences 
collected  en  June  7,  1982.  The  recorded  LORAN-C  rates  were 
also      9940-Y   and      9940-W.  Data   set      Number      2    is      located 

between  latitude  36o  40'  N  and  36°  45'  N  and  longitude  1210 
54'  M  and  122°  OO'W.  See  Figure  3.1  for  the  location  of 
test    areas    in   Monterey  Bay,    California. 

Data  set  Number  3  is  located  in  the  same  area  as  data 
set  Number  2,  between  latitude  36o  40'  N  and  36°  U5'  N  and 
longitude  1210  W  54»  and  122°  00'  W.  This  set,  consisting 
of  128  points  with  recorded  rates  9940-X  and  9940-Y,  was 
recorded  on  June   6,    1982. 

The  final  set,  data  set  Number  4  was  recorded  between 
June  3  and  iJune  5,  1982.  It  contains  169  points  located 
between  latitude  36o  36  N  and  36o  39'  N  and  longitude  121o 
53'  W  and  121°  58'  W.  The  recorded  rates  were  9940  -X  and 
9940-Y  which  are  the  same  as  data  set  three.  The  dana  for 
rates  9940-X  and  9940-Y  was  Icepr  separate  so  as  to  determine 
if  there  was  a  significant  difference  between  the  offshore 
(data  set  three)  and  inshore  (data  set  four)  drms  values  due 
to   phase   recovery    (see  Chapter   Two) . 
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TABLE    X 


Dara  Set   Parameter 


Set 
Number 

Collection 
Date 

TD 

R  ates 

N 
Da 

amber   of 
ta    Points 

Area 
Lit 

Limits 
Lon 

1 

June    12,    1979 
July    25,     1979 

9  940-Y 
9  9U0-W 

130 

36/38 
36/47 

121/49 
122/02 

2 

June   7,    1982 

9  940-Y 
9940-W 

193 

36/40 
36/45 

121/54 
122/00 

3 

June   6,    1982 

9  940-X 
9940-Y 

128 

36/40 
36/45 

121/54 
122/00 

4 

June   5,    1982 

9  94  0-X 
9940-Y 

169 

36/36 
36/39 

121/53 
121/58 
J 

Table  X  provides  a  convenient  breakdown  of  the  parame- 
ters for  each  of  the  data  sets.  The  table  consists  of  the 
data  set  number,  the  date  the  data  was  collected,  the 
LORAN-C  time  difference  rates,  the  number  of  data  poinxs, 
and  the  area  limits  in  latitude  and  longitude.  It  defines 
the  parameters  for  the  four  data  sets  of  Tables  XI  (Seawater 
Secondary  Factor  Errors)  ,  XII  (Semi-Empirical  TD  Grid 
Correction  Errors),  XIII  (Multiple  LOP.AN-C  Correc-ion  Table 
Errors) ,  and  XIV  (Multiple  Observed  Correction  Errors) .  All 
basic  data  and  computarion s  are  kept  on  file  with  the  NOAA 
hydrography  instructor  at  the  Naval  Pos-graduate  School, 
Monterey,  California.  (Individuals  seeking  this  information 
should  contact   the   Oceanography    Department.) 

A.       SEAWATER    SECONDARY    FACTORS 

LORAN-C  rates  were  computed  using  only  the  seawater 
Secondary  Factor  (equation  2. 13) .  The  LORAN-C  computed  time 
differences  are  basically  uncorrected  rates  since  no 
Additional    Secondary    Factors    were   applied.         Offsets    between 
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TABLE    XI 


Seawater  Secondary   Factor   Errors 


Rates 


Mean 
Offset 


Mean 
Offset 


Standard 

Error 
(s-^sec) 


a) 


b) 


c) 


d) 


Set 
Y 
W  . 

Set 
Y 
W 

Set 
X 
Y 

Se-c 
X 
Y 


(Schnebele's   data.    130    data    points) 
-0.508  -148.636  0.088 

-1.2U1  -683.693  0.134 

(June   7-    1982,    193   data    points) 
-0.526  -154.295  0.077 


-1.283 


-701.947 


0.  124 


(June    6,    1982,    128  data    points) 

1.565  535.449  0.059 


-0.550 


-161.404 


0.083 


(June    3-5,    1982,     169    data    points) 

1.582  550.089  0.055 

-0.817  -237.109  0.187 


drms 
(1) 


110.4 

101.3 

38.6 

68.  1 


the  observed  and  calculated  rates  using  only  the  seawater 
Secondary  Factor  were  computed  to  illustrate  -he  improvement 
in  posi-^,ion  after  applying  the  Semi-Empirical  TD  Model,  ASF 
Loran-C  Correction  Tables,  or  nhe  Multiple  Observed  ASF 
Correctors.  The  mean  offset  in  microseconds  and  meters, 
standard  error,  and  the  drms  are  found  in  Table  XI. 
Examples  of  offsets  for  several  data  points  are  listed  in 
Appendix   3    after   program   LORAN. 

It  was  stated  earlier  in  Chapter  One  that  Schnebele 
obtained  a  66  m  1  drms  using  Y  and  W  rates.  The  drms  of  66 
m  was  obtained  using  48  data  points  which  were  located  10  km 
or  more  offshore.  The  110.4  m  1  dris  for  the  Y  and  W  rates 
in  Table  XI  is  a  result  of  combining  the  130  inshore  and 
offshore  positions.  The  increase  from  66.0  m  1  drms  for  the 
offshore  positions  to  110.4  m  for  combined  offshore  and 
inshore    positions    indicatas    that   the   application    of   a    single 
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ASF  Corrector  is  dependant  on  the  size  of  the  area.  The  48 
data  points  were  located  between  latitude  36°  41 »  N  and  360 
46»    N   and   longitude    121o   55 •    W  and    122°    02«    W. 

Schnebele  also  obtained  the  66  m  1  drms  by  adjusting  the 
LORAN-C  observed  time  differences  which  were  skewed  due  to 
the  ship's  motion  and  the  five  second  averaging  interval  of 
the  LORAN-C  receiver.  These  caused  the  observed  TDs  to  be 
several  seconds  old  in  comparison  to  the  microwave  system 
measurements  [Schnebele,  1979].  Due  tc  the  large  amount.s  of 
data  from  the  June  1982  survey  operations,  no  deskewing  was 
done. 

The  offset,  standard  error,  and  drms  for  data  set  Number 
2  (June  7,  1982)  compares  well  to  data  set  Number  1 
(Schnebele's  data  -  all).  Also,  the  38.6  m  1  drms  for  the 
X-Y  rates  is  between  Schnebele's  U2.0  m  1  drms  prediction 
[Schnebele,  1979]  and  Nelson's  findings  of  30  m  1  drms  in 
San  Francisco  Bay  [General  Electric  3o. ,  1979].  The  large 
drms  of  68.1  m  for  data  set  Number  4  (X-Y  rates)  is  probably 
due  to  the  phase  recovery  of  the  electromagnetic  wave  from 
The  99aO-Y  transmitter  located  in  Searchlight,  Nevada. 
Finally,  if  ASF  Correctors  are  not  applied,  drms  values 
ranged  from  38.6  m  for  9940  X-Y  to  101.4  m  for  9940  Y-W  for 
data    sets      2   and   3.  The    large      offsets    for   all      data   sets 

indicate  a  systematic  error,  ranging  from  150  m  to  700  m, 
which  precludes  the  use  of  zhis  method  for  hydrographic 
surveying. 

B.       SEMI -EMPIRICAL    TD    GRID 

To  determine  if  the  Semi-Empirical  TD  Grid  would  reduce  the 
offset  and  drms  for  hydrographic  surveying,  program  LOPLC 
(Appendix  C)  was  applied  to  data  sets  1  through  4.  The 
Semi- Empirical  TD  Grid  applies  a  Secondary  Factor  and  an 
Additional    Secondary   Factor    to   the    primary    phase    delay    based 
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or.  the  distance  overland,  the  distance  over  water,  and  the 
total  distance  using  mean  land  and  water  conductivities. 
The  mean  offset  in  meters  and  microseconds,  standard  error, 
and  drms  are  listed  in  Table  XII.  Examples  of  offsets  for 
several  data  points  can  be  found  in  Appendix  C  after  program 

TABLE  XII 
Semi-Empirical  TD  Grid  Correction  Errors 


Rates 


Mean 

Offset 

(xz^sec) 


Wean 
Offset 
(111) 


Standard 
Error 


a) 


b) 


c) 


d) 


Set 
Y 
W 

Set 
Y 
H 

Set 
X 
Y 

Set 
X 
Y 


(Schnebele's   data,    130   data    points) 
1.131  330.672  0.086 


1.  145 


630.278 


(June   7,    1982  data,    193   data 

1.108  324.975 

1.105  603,926 

(June    6,    1982  data,    128    data 

1.131  386.940 


0.5  89 

(June  3-5, 
1.0  50 
0.341 


172.  923 


0.  139 

points) 
0.077 
0.  126 

goints) 
.061 
0.084 


1982    data,    169    data    points) 
365.113  0.052 

99.084  0.189 


drms 
(iS) 


113.5 


112.7 


39.4 


63.2 


LOPLC.  Earlier,  Program  LOPLC  had  been  tested  with  data 
found  in  the  article  by  The  Analytic  Science  Corporation 
[1979],  Program  LOPLC  produced  offsets  of  0.86  ;isec  and 
1.13  }is(^c  for  the  time  differences  from  rates  9940-X  and 
9940-Y,      respectively.  As    previously      stated,       the      large 

offsets  may  be  due  to  the  method  by  which  the  distance  and 
azimuth  were  computed.  The  drms,  for  the  data  from  Monterey 
Bay,  obtained  with  the  Semi -Empirical  TD  Grid  was  similar  to 
the  drms  errors  for  seawater  Secondary  Factor  .lodel.  This 
might    be   an  indication   that    Program   LOPLC    is   correct    but    the 
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bias  needs  to  be  adjusted  to  reduce  the  large  offset  between 
the      observed  and      calculated  TD      rates.  The   mean      offset 

ranged  from  99  to  630  m«  Bias  reduction  could  be  achieved  by 
applying  land-sea  data  for  the  entire  West  coast  to  the 
Semi-Empirical  model.  Again,  the  existence  of  large  offset 
values  precludes  the  use  of  this  method  for  hydrographic 
surveying. 

C.       TABLE    ASF    CORRECTORS 

To  determine  if  the  multiple  ASF  corrections  from  the 
LORAN-C  Correction  Table  would  reducs  the  offset  and  drms  to 
meex  the  NOS  accuracy  standards,  program  LORTAB  (Appendix  D) 
was    applied  to   data   sets    1      through   !4 .         ASF   Correctors   from 

TABLE   XIII 
Multiple   LORAN-C    Correction    Table   Errors 


Rates 


Mean 
Offset 
(iZ^sec) 


Mean 
Offset 

(iz^) 


standard 
Error 
(sz^sec) 


a)     Set    1     (Schnebele's   data,    130    lata    points) 


b) 


c) 


d) 


Y 
W 

Set 
Y 

Set 
X 
Y 

Set 
No 

eas 


-0.229 
0.  169 


-^7,068" 

92.540 


(June    7-    1982    data,    192   data 
-0.309  -90.614 


0.  123 


67. 189 


(June   6,    1982    data,    128   data 
0.559  19^.134 


-0.319 


-93.486 


0.  104 
0.  160 


oints) 
.085 
0,  144 


points) 
0.073 
0.  101 


arms 
(m) 


123.6 


1  16.9 


47.3 


4     (June    3-5,     1982    data,    169    data    points) 
ASF   Corrections   listed    in   Table    for   the    south- 
end    of    Monterey    Bay    next    to    the   shore    line. 


the      LORAN-C   Correction     Tables      are      determined    from      field 
observation  of   land   conductivities    by      the    U.S.      Coast    Guard 
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Calibration  Team  [Marine  Science  Department,  1982].  The 
results   are  listed   in    Tabls    XIII. 

When  compared  to  the  Saawater  Secondary  Factor  Error  in 
Table  XI,  ^he  drms  values  using  LORAN-C  Corrsction  Tables 
were  increased  slightly  while  the  offsets  were  reduced 
substantially.  The  offsets  ranged  from  65  -  200  m.  Since 
the  drms  ranged  from  47.3  to  116.9  m  for  rates  9940  X-Y  and 
9940  Y-W  respectively,  the  application  of  ASF  Corrector  from 
the  tables  does  not  meet  the  NOS  accuracy  standard.  Again, 
large  offsets  and  the  increase  in  drms  precludes  the  use  of 
this    method   for   hydrographic  surveying. 

D.       MULTIPLE    OBSERVED    CORRECTORS 

To  determine  if  multiple  observed  correctors  would 
diminish  the  offset  and  drms  values,  mean  ASF  Correctors 
were  selected  at  one  minute  latitude  and  longitude  in-ervals 
using  Program  ASFSEL.  One  minute  ASF  Correctors  are  shown 
at   the  end      of   Program  ASFSEL  in   Appendix   E.  From   the   one 

minute  grid,  mean  ASF  Correctors  were  selected  and  entered 
into  Subroutine  TABLE  at  five  minute  latitude  and  longitude 
intervals  in  Program  L0RTA3 .  The  following  is  an  example  of 
the  9940-X  ASF  Correctors  at  five  mi.iute  latitude  and  longi- 
tude   intervals    for   the  program: 

122/00/00.0  121/55/00.0  121/50/00.0 

36/50/00.0  -1.52  -1.63 

36/45/00.0  -1.52  -1.56  -1.61 

36/40/00.0  -1.58  -1.60 

This  is  the  same  format  used  in  the  LORAN-C  Correction 
Tables.  Program  LORTAB  was  applied  to  the  four  data  s*=ts. 
The    error   results   are   listed   in    Table    XIV. 

The  drms  value  obtained  with  multiple  observed  correc- 
tors   were      all   reduced      when   compared      to   the      irms    for      the 
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TABLE    XIV 
Multiple   Observed  Correct icn   Errors 


- 

Mean                    Mean               Standard 
Offset               Offset                 Error 

_   _   —     _      ^ 

drms 

Rates 

1 

(X-jisec)                l-ffl)                 (s-^sec) 

(Schnebels's   data,    130   data    points) 
0.033                    9.585                   5.076 
0.052                28.604                   0.113 

(I) 

87.6 

a) 

Set 
Y 
W 

b) 

Set 
Y 

2 

(June   7,    1982,    193  data    poinds 

0.028                    8.318                   0.073 
-0-016                -8. 944                   0.116 

89.3 

c) 

Set 
X 
Y 

3 

(June   6,    1982,    128  data    points) 
-0.034                -11.490                 0.055 
0.010                     2.913                0.086 

38.3 

d) 

Set 
X 
Y 

a 

(June    3.5,    1982,     169    data    points), 

0.001                     0.383                0.052 

-0.008                   -2.225                 0.187 

67.5 

seawater  Secondary  Corrector.  The  idst:  impressive  reduction 
in  drms  was  within  Schnebele's  data  which  covered  an  area  of 
seven  minutes  of  latitude  and  12  minutes  of  longi-ude.  The 
drms  for  seawater  Secondary  Correctors  was  110.4  m  whereas 
the  drms  for  the  same  data  using  multiple  observed  correc- 
tors was  87.6  m.  This  is  a  smaller  drms  than  that  of  the 
June  7,  1932  data  (data  set  Number  2)  which  was  obtained 
three  years  later.  It  appears  that  87.6  m  1  drms  is  nearly 
the  minimum  error  that  can  be  obtained  for  the  9940  Y-W 
rates  in  Monterey  Bay  after  applying  multiple  observed  AS? 
Correctors  at  five  minute  latitude  and  longitude  intervals. 
For  the  LORAN-C  rates  9940  X-Y,  a  38,  3  m  1  drms  was  obtained 
for  same  five  minute  area  covered  by  data  set  Number  2  (June 
7,  1982).  LORAN-C  rates  9940  X-Y  were  not  obtained  for  the 
same  size  area  covered  by  data  set  Number  1  (seven  minutes 
of   latitude   and    twelve   minutes    of   longitude   -    9940      Y-W)     due 
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to  the  restriction  of  ship  time  and  the  length  of  time  the 
Racal-DECCA  Trisponder  alrctronic  eqaipinent  had  been  loaned. 
The  drms  for  the  June  3-5,  1982  inshore  data  was  only 
reduced  to  67.5  m  from  68.1  m  for  seawater  Secondary 
Correctors.  The   small      change    in      error   at      the   coast      is 

probably  a  result  of  the  erratic  behavior  of  phase  recovery 
discussed  earlier  in  Chapter  Two.  The  drms  value  could 
presumably  be  reduced  if  the  correctors  were  applied  at  one 
minute   intervals.  This   would    be     a    very   costly      method   of 

calibrating  Loran-C   for   hydrographic   surveying. 
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V.       CONCLUSIONS 

It  was  noted  in  Chapter  One  that  the  smallest  scale 
routinely  used  for  coastal  surveys  Ls  1  ;80,000.  This  yields 
an  allowable  error  of  49  tn  1  drms  with  no  systematic  errors. 
This  paper  determined  whether  or  not  applying  multiple 
Additional  Secondary  Factors  (ASF)  Correctors  to  LORAN-C 
lines  of  position  would  reduce  the  drms  sufficiently  to  meet 
the   accuracy   standards  set    by  the   National   Ocean    Survey. 

Three  methods  of  applying  multiple  ASF  Correctors  were 
tested.  The  first  approach  computes  the  time  difference 
based  on  a  Semi-Empirical  ID  Grid.  The  Semi-Empirical  Xodel 
produced  large  offsets  in  the  9940-W,  99aO-X,  and  99aO-Y 
time  differences.  The  offsets  ranged  from  99  to  630  m.  The 
drms  for  9940  X-Y  combination  was  39.4  m  and  the  drms  for 
9940    W-Y   combination    was    102.7  m. 

The  second  method  applies  AS?  Correctors  found  in  the 
DMAHTC  LORAN-C  Correction  Tables  to  LORAN-C  lines  of 
position.  The  application  of  the  tables  reduced  the  offset 
in   the   LCRAN-C   time  differences.  The   offsets   were    between 

67  and  191  m.  The  drms  was  increased  to  47.3  m  for  the  9940 
X-Y    combination    and    116.9    m    for    the    9940    W-Y    pair. 

The  final  and  most  accurate  method  applies  multiple 
observed  ASF  Correctors  at  five  minute  latitude  and 
longitude      intervals    to      LORAN-C   lines      of    position.  This 

method  again  reduced  the  offset  ia  the  time  difference. 
This  offset  was  between  3  and  12  m  for  the  9940  X-Y 
combination.  Part  of  the  offset  may  have  been  a  result  of 
the  microwave  positioning  system.  Reference  is  made  to 
Anderson's  [1982]  paper  (in  preparation)  on  the  evaluation  of 
multiple   lines   of    position. 
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The  drms  values  were  also  reduced  to  38.3  m  for  the  99U0 
X-Y  rates  and  89.3  m  for  the  9940  W-Y  combination.  The  38.3 
m  1  drms  can  be  decreased  by  improving  the  sampling  time  for 
LORAN-C  receivers.  Nelson  obtained  30  m  1  drms  for  the  9940 
X-I  rates  with  special  LORAN-C  squipment  used  in  San 
Francisco  Bay,  California  [General  Electric  Co.,  1979]. 
Improving  the  sampling  time  for  LORAN-C  receivers  used  as 
positioning  equipment  for  hydrographic  surveys  should  be 
investigated. 

With  drms  values  of  38.3  m  with  the  possibility  of 
obtaining  30.0  m  1  drms  and  offsets  ranging  from  3  to  12  m, 
it  may  be  possible  to  use  LORAN-C  for  hydrographic  surveys 
at  scales  of  1:80,000  or  less  using  multiple  observed  ASF 
Correctors.  The  use  of  DMAHTC  LORAN-C  Correction  Tables 
should  not  be  ignored.  After  updating  these  ASF  Correctors 
with  observed  data,  the  LORAN-C  Correction  Tables  may  allow 
LORAN-C  to  be  used  as  a  positioning  system  for  hydrographic 
surveys. 
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APPENDIX    A 
PROGRAM    GPBYLQ 

C  PROGRAM    GPBYLQ 

C 

C  GENERAL    PROGRAM    FOR    DETERMINING    GP    FROM    KNOWN    STATION 

C  POSITIONS    AND    THE    DISTANCES    FROM    THEM    USING    LEAST 

C  SQUARES.       GP    TO    XY    AND    XY    TO    GP    ARE    DETERMINED    3Y 

C  SUBROUTINES    GPTOXY    AND    XYTOG?    WHICH    ARE    BASED    ON    THE 

C  MODIFIED    TRANSVERSE    MSRCATOR    PROJECTION     (MTM) . 

C 

C  PROGRAMMED  BY  GERALD  E.  WHEAION,  LT.  NOAA 

C 

C  LEAST  SQUARES  ADJUSTMENT  BY  PAUL  R.  WOLF,  PH.D.  AND 

C  REVISED  BY  D.  LEATH,  LCDR 

C 

C  PROGRAM    INPUT    VARIABLE    NAMES 

C  TITLE    =    ANY    JOB    IDENTIFICATION    NAMES    OR    NUMBERS 

C  M    AND    N    =    THE    NUMBER     OF    EQUATIONS     (M)     AND    UNKNOWNS     (N) 

C  XO    AND    YO    =    BEST   QUESTIMATE    OF    THE    POSTION 

C  STA(I,1)     =    X    COORDINATE    OF    KNOW    STATION 

C  STA(I,2)     =    Y    COORDINATE    OF    KNOW    STATION 

C  STA(I,3)     =    DISTANCE    FROM    KNOW    STATION 

C  A  (I,  J)     =    THE    COEFFICIENT    MATRIX 

C  EL  (I,  J)     =    THE    CONSTANT    MATRIX 

C  QLL{I,J)     =    THE    WEIGHT    MATRIX     (WEIGHTS    ARE    ENTERED    AS 

C  1 • S    IF    THE    SOLUTION    IS    EQUALLY    WEIGHTED) 


C 


IMPLICIT    REAL*8     (A-H,0-Z) 

DIMENSION    XCORD  (30) ,Y CORD (30)  ,WT  (3  0)  ,ISNO(30) 

COMMON    /ISTAT/    STA  (30  ,  3)  ,  FXY  (30) 

COMMON    /LSQX/    A  (30  ,30  )  ,  EL  ( 30  ,  1)  ,  QLL  ( 30  ,  30)  ,  AT  (30  ,  30)  , 
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1AQ(30,30),    QXX(30,30)   ,    AQL(30,10),    X(30,1),    V(30,1), 

2VAE(30)  ,    TITLE(80) 
C 

C  READ    AND    WRITE    OaTPOT    TITLE 

C 

WRITS  (6,509) 

READ  (5,710)  TITLE 

WRITE  (6,710)  TITLE 
C 

C  DEFINE    NUMBER    OF    SI3NALS     (NOT    GREATER    THAN    30)     AND 

C  NUMBER    OF    DATA    SETS. 

C 

READ  (5,502)     NSIG,NDATA 

502  FORMAT(I5,I10) 
C 

C  DEFINE    THE    CENTRAL    MERIDIAN     (CMSR) ,     FALSE    EASTING 

C  (FEST),    AND    CENTRAL    LATITUDE     (CLAT)  . 

C  CMER    IS    EXPRESSED    IN     DEGREES,     MINUTES,    AND    SECONDS. 

C'  FEST    IS    THE    X-COORDINATS    THAI    IS    ASSIGNED    TO    THE 

C  CENTRAL    MERIDIAN    AND    IS    EXPRESSED    IN    METERS. 

C  CLAT    IS    DEFINED    AS    THE    CONTROLLING    LATITUDE. 

C  IT    IS    USED    TO    REFERENCE   THE    Y-COORDIN ATES    AND 

C  IS    EXPRESSED    IN    METERS. 

C 

READ  (5,503)  ILO NC, ILM INC, RLSECC , FEST , CLAT 

503  F0eMAT(1X,I3, 1X,I2,1X,F8. 5, 1X,?7. 1,F10.1) 

CMER    =     ((IA3S  (ILONC)      *    60     +    ILMINC)     *    60)     +    RLSECC 
C 

C  DEFINE    VARIABLE    FOR: 

C  NUMBER    OF    EQUATIONS    (M) 

C  NUMBER    OF    UNKNOWNS     (N) 

C  IPAGE    =    NUMBER    OF    LINES    PER    PAGE. 


C 


IPAGE    =    1 


66 


M  =  4 

N  =  2 
C 

C      NULL  WEIGHTS 
C 

DO  3   I a= 1,30,1 

DO  3   JM=1,30,1 
3      QLL(IM,JM)  =0.0 
C 

C      R2AD  STATION  NUMBERS  , POSITION,  AND  WEIGHTS. 
C      CONVERT  POSITIONS  (3P)  TO  SECONDS  AND  THEN  TO  XY. 
C 

DO  12   J=1,NSIG,1 

READ  (5,80  0)  I SNO  (J)  ,IL  AT  ,1  MIN,  RS  EC,  JLON,  JMIN  ,  SSEC,  WT  (J) 
800         FORMAT(1X,I3,IU,I3,F9 .5,I5,I3,F9.5,F5. 1) 
C 

RMAST    =     ((lABS(ILAT)     *    60    +    IilIN)     *    60)     +    RSEC 

RMASTL    =     ((lABS(JLON)      *    60    +    JMIN)     *    60)     +    SSEC 
C 

CALL    GPTOXY (RMA ST,  RMA STL, XMETER, YMETER, FEST , CLAT , CMER) 

XCORD(J)     =    XMETER 
12  YCORD(J)     =    YMETER 

C 

C  READ    DATA        (STATION    NUMBERS    AND    THE    DISTANCES) 

C 

15  DO    40      JC0UNT=1  ,NDATA  ,1 
READ(5,805)     I F, IS, IT , lU ,NREC 

805  FORM  AT  (515) 

DO    16       IC0UNT=1  ,NSI3,  1 
IF  (ISNO  (ICOUNT)     .EQ.     IF)     IF=ICOUNT 
IF  (ISNO  (ICOONT)     .EQ.     IS)     IS  =  ICOUNT 
IF  (ISNO  (ICOUNT)     ,EQ.     IT)     IT=ICOUNT 

16  I?  (ISNO  (ICOUNT)     .EQ.     14)     I4  =  IC0UNT 
C 
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STA(1,1)     =    XCORD(IF) 

STA(1,2)     =    YCORD(IF) 

QLL(1,1)     =    WT(IF) 
C 

STA(2,1)     =    XCOFD(IS) 

STA(2,2)    =    YCORD(IS) 

QLL(2,2)     =    WT(IS) 
C 

STA(3,1)     =    XCORD(IT) 

STA(3,2)     =    YCORD(IT) 

QLL(3,3)     =    WT(IT) 
C 

STA(a,1)     =    XCORD(IU) 

STA(4,2)     =    YCORD(IU) 

QLL(4,4)     =    WT(I4) 
C 

C  READ    THE    DISTANCE    RECORD    AND    LORAN    HATS. 

C 

DO    38    KC00NT=1, NREC,1 

READ  (5,507)     STA  (1,  3)  ,  STA(2,3)  ,STA  (3,3)  ,STA  (4,3)  , 
1    RATE1,RATE2 
507  FORMAT(6F10,  1) 

C 

C  DETERMINE    BEST    GUESS     COORDINATES    XO    AND    YO 

C  WITH    SUBROUTINE    GUESS. 

C 

19  CALL    GUESS (XO,YC) 
C 

IJUMP    =    0 
C 

C  COMPUTE    FXY,     A    AND    L     MATRIX 

C 

20  DO    25    1=1, M 

FXY(I)     =    DSQRT  (DA3S(  (  XO-STA  (I,  1)  )  **2 
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1  (Y0-STA(I,2)  )**2)  ) 

A  (1,1)  =  (X0-STA(I,1)  )  /  (FXY(I)) 

A(I,2)  =  (Y0-STA(I,2)  )  /  (FXY  (I)  ) 
25     EL  (1,1)  =  STA(I,3)  -  FXY(I) 
C 

C      CALL  SUBROUTINE  LSQR 
C 

CALL  LSQR  (K,N) 
C 

C      COMPUTE  THE  NEW  QUESSTIMATE  FOR  XO  AND  YO 
C 

XO  =  XO  +  X  (1  ,1) 

YO  =  YO  +  X  (2,1) 
C 

C      EXIT  IF  STANDARD  ARE  MET  USING  IJUMP  OR 
C      XO  AND  YO  CUT  OFF 
C 

IF  (DABS  (X  (1,1)  )  .LE.1.  00  .AND.  D  AB3  (X  (2,  1)  )  .  LE.  1  .  00) 
1   GO  TO  35 

IF  (IJUMP  .EQ.  10)  GO  TO  35 
30     GO  TO  20 
C 

C      COMPUTE  ERROR  ELLIPSE 
C 
35     CONTINUE 

CALL  ELIPSE  (SU,  SV) 
C 

C  CONVERT    XY    TO    GP 

C 

CALL    XYTOGP (XO, YO, SEC lAT, SECLON , FEST , CLAT, CMER) 

CALL    TODMS (SECL AT, IDE GP,I MINP, RSECP) 

CALL    TODMS (SECLON, JDE GP ,JMINP , 3SSC?) 
C 
C  PAGE    AND    CONTINUE    WITH    NEXT    SET    OF    OBSERVATIONS 
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c 

WRITE  (6,505)     IDEGP ,IM INP,RSEC?, JDEGP, JMINP , SSECP , 
1    EATE1,RATE2 
C  WRITE(6,505)     IDEGP ,IMINP, RSECP, JDEGP, JMINP, SSECP, 

1  sa,sv 

505         FORMAT (14, 13, F6. 2, 15, I3,F6 . 2, 2F9 . 2) 
C  IF(IPAGE    .EQ.    50)     WRITE  (6,509) 

IF(IPAGE    .EQ.     50)     IPAGE   =    0 
38  IPAGE    =    IPAGE    +    1 

C 

C      FORMAT  STATEMENTS 
C 

710         FORMAT(80A1) 
509  FORMAT  (1H1) 

UO  CONTINUE 

STOP 

END 

C========  =  =======  =  ==  ======  ======  =  ==  =  =  =  ==  ======  =  =  ==  =  ======: 

SUBROUTINE    LSQR  (M,N) 
Q 

IMPLICIT    REAL*8  (A-H,0-Z) 

COMMON    ASQX/    A  (30,30)  ,SL  (30,  1)  ,QLL  (30,30)  ,  AT(30,30)  , 
1AQ(30,30),    QXX(30,30)  ,    AQL(30,10),    X(30,1),    7(30,1), 
2VAR(30)  ,    TITLE  (80) 
C 

C  COMPUTE    A    TRANSPOSE    3Y    TRANSPOSING    THE    A    MATRIX     (AT) 

C 

DO    6  1    I=1,M 
DO    61    J=1,N 
61    AT  (J, I)  =A  (I, J) 
C 

C      USING  STEPS  (1)  ,  (2)  ,  AND  (3)  COMPUTE  THE  INVERSE 
C      OF  THE  TRANSPOSE  (AT)  *  WEIGHTED  MATRIC  (QLL)  * 
C      MATRIX  A  =  QXX. 
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c 

C  (1)        COMPOTE    AQ   =    AT    *    QLL 

C 

DO    7  1    1  =  1 ,N 

DO    71    J=1 ,M 

AQ  (I,J)=0. 

DO    71    K  =  1  ,M 
71    AQ  (I,J)=AQ(I,J)  +(AT(I,K)*QLL(K,  J)  ) 
C 

C  (2)        COMPUTE    QXX    =    AQ    *    A 

C 

DO    81    1  =  1,  N 

DO    81    J=1,N 

QXX(I,J)=0. 

DO    81    K=1,M 
81    QXX(I,J)=QXX  (I,  J)  +AQ  (I,K)  *A  (K,J) 
C 
C  (3)        INVEST    QXX    MATRIX 


C 


DO    3  07    K=1 ,N 
DO    302    J=1,N 
IF  (J-K)  304,302,  304 
304    QXX(K,J)=QXX  (K,  J)/QXX  (K,K) 

302  CONTINUE 
QXX(K,K)=1./QXX  (K,K) 
DO    307    1=1, N 

IF  (I-K)  305,307,305 
3  05    DO    3  03    J=1,N 

IF  (J-K)  306,303,  306 

306  QXX(I,J)=QXX  (I,  J)  -QXX  (I,K)  *QXX  (K,  J) 

303  CONTINUE 
QXX(I,K)=-QXX  (I,K)  *QXX  (K,  K) 

307  CONTINUE 
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C      USING  STEPS  (4)  AND  (5),  COMPOTE  THE  UNKNOWNS  X 

C      BY  MULT  THE  INVERSE  Q XX  AND  AQL. 

C 

C       (U)   COMPUTE  AQL  =  A  Q  *  EL 

C 

DO  101  1=1,  N 

AQL(I,1)=0. 

DO    101    K=  1 ,  M 
101    AQL(I,1)=AQL  (I,  1)  +A3  (I,K)  *EL(K,  1) 
C 
C  (5)        COMPUTE    X    =    QXX    *    AQL 

DO    201    1=1,  N 
C 

X  (I,  1)=0. 

DO    201    K=1,N 
201    X  (I,  1)  =X(I,1)  +QXX(I,K)  *AQL  (K,  1) 
C 

C       (6)   COMPUTE  THE  RESIDUAL   (V  =  A  *  X  -EL) 
C 

DO  301  1=1,  M 

V{I,1)=0. 

DO  3  01  K=1,N 
301  7  (I,  1)=7(I,  1)  +A  (I,K)  *X(K,  1) 

DO  1  1=1, M 
1  V(I,1)=V(I,1)-EL(I,1) 
C 

C      COMPUTE  THE  STANDARD  DEVIATION  OF  UNIT  WEIGHT  SIGMA 
C         DM  -  NUMBER  OF  OBSERVATIONS 
C         DN  -  NUMBER  OF  KNKNOWNS 


C 


SIGMA=0. 

DM=M 

DN=N 

DO  332  1=1, M 


72 


382    SIGMA=SIGMA+V  (1,1)  **2*QLL  (1,1) 
SIGMA=DSQRT     (SIGMA/(D M-DM) ) 
C 

C  COMPUTE   THE    STANDARD     DEVIATION    OF    THE    ADJUSTED    UNKNOWNS 

C  QXX    -    ARE    THE    ELEMENTS    OF    THE    COVARIANCE    MATRIX. 

C 

DO    446    1=1, N 
446    VAR(I)  =DSQRT  (QXX(I  ,1)  *SIGMA**2) 
C 

5  10    CONTINUE 
RETURN 
END 

C==r==  =====  =  ==  ==============  =  =  ========  =  =  =  =  ==  =  =  =  ===  =  ==  =  =  ======  = 

SUBROUTINE    ELIP SE (SU, SV) 
C 

c 

C  SOLVE    FOR    THE    SEMI    MAJOR    AND    SEMI    MINOR    AXIS    OF 

C  THE    ERROR    ELLIPSE 

C 

IMPLICIT    REAL*8  (A-H,0-Z) 

COMMON    /LSQX/    A  (30  ,30  )  ,  EL  (  30,  1)  ,  QLL  (30  ,  30)  ,  AT  (30  ,  30)  , 
1AQ(30,30),    QXX(30,30)  ,    AQL(30,10),     X(30,1),     V(30,1), 
2VAR(30)  ,    TITLE(80) 
C 

SUS    =    .5*(QXX(1,1)     +     QXX(2,2)     +    DSQRT  (DABS  (QXX  ( 1 ,  1)  - 
1    QXX  (2,2)  +4.0*QXX(1,2)  *QXX  (2,  1)  )  )  ) 

SVS    =    .5*(QXX(1,1)     +     QXX(2,2)     -    DSQRT  (DABS  (QXX  (1  ,  1) - 
1    QXX  (2,2)  +4.0*QXX(1,2)*QXX  (2,1)  )  )  ) 
C 

SU    =      DSQRT  (SUS) 
SV    =       DSQRT  (SVS) 
C 

RETURN 
END 
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C==:==s  =  =:  ==  =  =====  =  =  =  =  =======  =====  =  =  =  =  =  =  =  ===  =  ==  =  ==  =  =  ==  =  =  =  =  ==  = 

SD3R0UTINE    G0ESS(XO,YO) 

c 

C  SUBROUTINE    GUESS   DETERMINES    THE    BEST    GUESS    COORDINATES 

C  TO    BE    USED    IN    SUBROUTINE    LSQR.       USE    RIGHT    SIDE    RULE    FOR 

C  STATION    ORDER. 

C 

IMPLICIT    REAL*8     (A-H,0-Z) 

COMMON    /ISTAT/    STA  (30  ,  3)  ,  FXY  (30) 
C 

C  DETERMINE    DISTANCE    BETWEEN    STATION    1    AND    STATION    2 

C 

D  =  DSQRT  ((STA  (2,  1)  -STA  (1,1)  )  **2+  (STA  (2,2)-3TA  (1,2)  )**2) 
C 

C  DETERMINE    ANGLE    ALPHA    BETWEEN    XO , Y0/STA2/STA1 

C 

ALPHA=DARC0S ( (STA(2,3) **2-STA (1,3)  **2  +  D**2) / 
1     (2.0*STA  (2,3)*D)  ) 
C 

C  DETERMINE    ANGLE    BROVO    BETWEEN    X-AXIS    AND    STA2-STA1 

C 

BROVO=DARSIN  (  (STA  (1,2)-STA  (2,2)  )  /D) 
C 

C  DETERMINE    X    AND    Y    LENGTH 

C 

X=STA  (2,3)  *DCOS  (ALPHA +BROVO) 

Y  =  STA  (2,3)  *DSIN  (ALPHA  +BROVO) 
C 

C  DETERMINE    XO    AND   YO 

C 

XO    =    STA(2,  1)  +X 

YO   =    STA(2,2)  ■^Y 


in 


RETURN 
END 


SUBROUTINE    GPTOXY(SEC LAT, SECLON,XCO, YCO^FEST ,CLAT,CMER) 


IMPLICIT    REAL*8  (A-H,0-Z) 

DATA    E2,RK0,A    / .00676 8658  DO ,.  99998D0 , 6378206 .UDO/ 

DATA    RKGE0,W1,W2   /O. 0 4848 1 368D0 , 0. 11 422D0, 2 1 . 73607D0/ 

DATA    W3, W 4/5 104. 57 33 8  DO, 63 673 99. 6 3 9D0/ 

DATA    RADSEC    /.O 000048 48  13681 1 1D0/ 


C 
C 


RADLAT    =    SECLAT    *    RADSEC 
SINLAT    =    DSIN  (RADLAT) 
SIN2LA    =    SINLAT    *    SINLAT 
COSLAT    =    DCOS  (RADLAT) 
C0S2LA    =    COSLAT    *    COSLAT 
P    =     (CMER    -    SECLON)     /     10000. ODO 
V    =    A    /   DSQRT(1.0D0    -    E2*SIN2LA) 
TANCON    =    1.0D0    -    SIN2 LA/C0S2LA 

S    =    W4    *     (RADLAT    -    SI NLAT*C0SLAr/1 0. 0DO**6    * 
*  (W3    -COS2LA*(W2- W1*C0S2LA)  )  ) 

T1    =    S    *    RKO    -    CLAT 
T2    =    RKGEO    *    COSLAT    *    RKO    *    V 
T3   =    T2    *    RKGEO    /   2. 0  DO 
T4    =    T3    *    SINLAT 

T5    =    T3    *    RKGEO    *    COS 2LA   /    3. ODO 
T6    =    T5    *    TANCON 

T7    =     (4. ODO    +    TANCON)      *T5    *RKGEO    *    SINLAT    /    4, ODO 
XCO    =     (T2    +     (T6*P**2))*P    +    ?EST 
YCO    =     (T7*p**4)     +     (T4*P**2)     >    T1 
RETURN 
END 
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SUBROUTINE    XYTOGP(XCO , YCO, SECLAT, SECLON,FEST,CLAT,CMER) 


IMPLICIT    REAL*8     (A-H,0-Z) 
C 

DATA    E2,A,SR    /O .00676 86 58D0 ,6378206. UDO , 
1    0.0000048U81368D0/ 

DATA    W1,W2,W3    /O. 2468 2D0, 30. 0233 5D0, 5078 .64 97700/ 
C 

D    =    CLAT    ♦    YCO 

WO    =    0.  15704998  1D0/10.0D0**6    *    D 

SINWO    =    DSIN(WO) 

COSWO    =    DCOS  (WO) 

C0S2W0    =    COSWO    *   COSWO 

PHI1     =    WO    +    SINWO*COS WO/1 0.0D0**6     * 
*  (W3-^COS2WO*{W2+W  1*C0S2W0)  ) 

PHI2    =    PHI1    /   0.99998D0 

PHI3    =    PHI2    /   SR 

Q    =     (XCO    -    FEST)     /    10.0D0**6 

V    =    A    /   DSQRT{1.0D0    -    E2*DSIN  {PHI2)  **2) 

T    =    DC0S(PHI2)     *    SR 

C    =    V    *    0.99998D0 

T1    =     10.0D0**6    /     (T*C) 

T2    =     (T1*10.  0DO**6)     /     (2.0D0*C) 

T3    =     (T2*10.0D0**6)     /     (3.0D0*C) 

T4    =     (T3*10.0D0**6)     /     (4.0D0*C) 

DELLON    =     (T1-Q**2*T3*  (2.0D0*DTAN (PHI2) **2+1 .ODO) ) *Q 

SECLAT    =     (  (  (3.0D0*DTAN(PHI2)  **2  +  5.0D0)  *T4)  *Q**2-T2) 
1         *Q**2*DSIN  (PHI2)     +    PHI3 

SECLON    =    CMER    -    DELLON 

RETURN 

END 
SENTRY 
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DATA    SET    EXAMPLE      -       PROGRAM    GPBYLQ. 

SHIP'S    POSITION  OBSERVED    LORAN    RATES 


LATITUDE 
(D-M-S) 

LONGITUDE 
(D-M-S) 

9940-X 
(fisec) 

9940-Y 

3  6    36    40.20 

121    52    48,62 

27508.79 

42742.71 

36    36    42.41 

121    52    48.16 

27508.81 

42742.93 

3  6    36    44.56 

121    52    47.52 

27508.83 

42743.15 

3  6    36    46.84 

121    52    46.98 

27508.76 

42743. 26 

36    36    49.09 

121    52    46.5  9 

27508.76 

42743.48 
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APPENDIX    B 
PROGRAM    LOR  AN 

C  PROGRAM    LORAN 

C 

C  PROGRAM    COMPOTES    LINE    OF    POSITIONS    FOR    LORAN-C    USING 

C  THE    SECONDARY    FACTOR     (SF)     BASED    ON    SEAWATER    EM    MODEL. 

C 

IMPLICIT    REAL*8     (A-H,0-Z) 

DIMENSION    AXIS(  13)  ,RF  (13) 

DIMENSION    XSLAV1 (100  0) ,XSLAV2 (1000) 

DATA    RHOSEC,PI,UNCOV, RN/2. 06264 8062U71D05, 
1    3. 1 4 15 92 6535 89 8D0, 2 9  9.79  245  800,1 .000338DO/ 

DATA    XMEAN1,XMEAN2,VAR1,7AR2/O.OODO,O.OODO, 
1    0.0ODO,00.0D0/ 

DATA    XMEAN3,XMEAN4/0. 0ODO,O.00DO/ 
C 
CCCCCC* ********************  ************* 

DATA    AXIS/6. 3782064D06, 6. 378388         D06 , 6. 377397 1 55D06  , 

1  6.37816    D06, 6. 37816D06, 6, 378249145006,6. 378165D06, 

2  6.378166006,6.378165    D06 , 6 . 378 1 45D06 , 6. 3775634D06, 

3  6.378245006,6.3781350006/ 

DATA    RF/6. 3565838  006,2.97         002,2.9915  28  12  85002, 

1  2.98  2500  2,2. 98  2  4716  7  4270  02,2.9  3  46  50  02,2.9  8  25002, 

2  2.983    002,2.983002,2-9825    002,6.3562569006, 

3  2.933    002,2.9826002/ 
C*****ELLIPSOID    OPTION    NUMBER 

C  1.  CLARKE    1866  8.  MERCURY 

C  2.  INTERNATIONAL    (HAYFORD)  9.  MARSHALL    ISLAND 

C  3.  BESSEL     1841  10.  NAVY    8D 

C  4.  AND     (AUSTRALIAN)  11.  AIRY 

C  5,  1967    REFERENCE  12.  KRASSOWSKI     1940 
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C  6.    CLARKE    1880    MOD  13.     WGS    1972 

C  7.    SAO 

C 

cccccc* ************ ********  ************* 

C  CC1-2    =    ELLIPSOID    NUMBER    (K) 

C  CC3-5    =    NUMBER    OF    POINTS    ALONG    COAST     (IREC) 

C 

READ  (5,  100)     K,IREC 
100  F0RMAT(lX,l2,ia) 

TW0PI=2.*PI 

A  =  AXIS(K) 

F=1./RF  (K) 

IF  (F.LT.3.D-3)     F=(A-1./F)/A 
C 

C  READ    MASTER    AND    SLAVE    STATIONS    POSITIONS 

C  THE    FIRST    RECORD    IS    THE    NUMBER    OF    MASTER    AND    SLAVE 

C  STATIONS    FOR    THE   PARTICULAR    CHAIN. 

C 

WRITE  (6,202) 
202  FORMAT(IHI) 

WRITE  (6,201) 
201  FORMAT  (1H    ) 

READ  (5, 105)     ILATM,IMINM,RSECM,ILONM,ILMINM,RLSECM 

WRITE (6,1 05)     ILATM,IMINM,RSECM,ILONM,ILMINM,RLSECM 
105         FORMAT (IX, 13, IX, 12, IX, F5. 2, IX, 14, IX, 12, IX, F5. 2) 

RMAST    =     (  (lABS(ILATM)      *    60    +    IMINM)     *    60    >    RSECM)     / 
1    RHOSEC 

IF(ILATM    .LT.    0)     RMAST    =    -RMAST 

RMASTL    =     (  (TABS  (ILONM)     *    60    +    ILMINM)     *60    +RLSECM)     / 
1    RHOSEC 

IF  (ILONM    .GT.    0)     RMASTL   =    TWOPI    -    RMASTL 
C 

C  READ    THE    FIRST    SLAVE    STATION    AND    CHANGE    THE    LATITUDE 

C  AND    LONGITUDE    INTO    RADIANS. 
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READ  (5, 101)     ILATS,IMINS,RSECS,rLONS,ILMINS,ELSSCS, 
1    DELAY1 

WRITE  (6,101)     ILATS,IMINS,RSECSrILONS,ILMINS,RLSBCS, 
1    DELAY1 
101         FORMAT (IX, 13, IX ,12 , 1 X ,F5. 2, IX, lU , IX, 12, IX,  F5 . 2 , F9 . 2) 
RSLAV1    =     ((lABS  (ILATS)     *    60    4-    IMINS)     *    60    +    RSECS)     / 
1    RHOSEC 
IF  (ILATS    .LT.    0)     RSLAV1    =    -RSLAV1 

RSLAL1    =     ((lABS  (ILONS)     *    60    +    ILMINS)     *    60    ■»•    RLSECS)     / 
1    RHOSEC 
IF  (ILONS    .GT-    0)     RSLAL1    =    TWOPI    -    RSLAL1 
C 

C  READ    THE    SECOND    SLAVE    STATION    AND    CHANGE   THE    LATITUDE 

C  AND      LONGITUDE    INTO    RADIANS. 

C 

READ  (5,  101)     ILATS,  IMINS,  RSECS,  ILONS,  ILMINS,  RLSECS, 
1    DELAY2 

WRITE  (6,101)     ILATS, IMINS, RSECS, ILONS, ILMINS, RLSECS, 
1    DELAY2 

RSLAV2    =     ((TABS  (ILATS)     *    60    ^    IMINS)     *    60    ■••    RSECS)     / 
1    RHOSEC 
IF  (ILATS    .LT.    0)     RSLA V2    =    -RSLA72 

RSLAL2    =     (  (IA5S  (ILONS)     *    60    +    ILMINS)     *    60    +    RLSECS)     / 
1    RHOSEC 
IF  (ILONS    .GT.     0)     RSLAL2    =    TWOPI    -    RSLAL2 
WRITE  (6,201) 
C 

C  READ    THE    RECORDS   POSITIONS    AND    OBSERVED    RATES,     THEN 

C  COMPUTE    THE    TOTAL    DISTANCE. 


C 


DO    550    I=1,IREC,1 

READ (5, 104)     ILAT,IMIN ,RSEC ,ILON , ILMIN, RLSEC, 
1    RATE1,RATS2 


30 


104         FOEMAT(1X,I3,I3,F6.3,  15,13,  F6  .  3  ,  3X,  2F1  0  .2) 

C 

C  CONVERT    LAT    AND    LONG    TO   RADIANS. 

C 

RPOST    =     ((lABS(ILAT)      *    60    +    IMIN)     *    60    +    RSEC)     / 
1    RHOSEC 

IF(ILAT    .LT.     0)     RPOST    =   -RPOST 

RPOSTL    =     ((lABS(ILON)      *   60    +    ILMIN)     *    60    +    RLSEC)     / 
1    RHOSEC 

IF(ILON    .GT.     0)     RPOSTL    =    TWOPI    -    RPOSTL 
C 

C  COMPUTE    DISTANCES    AND    AZIMUTHS    FROM    THE    OBSERVED    POINT 

C 

C  MASTER 

C 

CALL    INVER1 (A, F, RPOST , RPOSTL, RM AS T, RMASTL, FAZM , BAZM , 
1    DISTM) 

DDISTM    =     (RN    *    DISTM)     /    UNCOV 

CALL    SECFAC  {UDISTM,SFM) 
C 

C  SLAVE1 

C 

CALL    INVER1 (A,F ,RPOST , RPO STL,RSLA V1 , RSLAL1 , F AZ 1 , SAZ 1 , 
1    DIST1) 

UDIST1    =     (RN    *    DIST1)     /   UNCOV 

CALL    SECFAC  (UDI  ST1, SF  1) 
C 
C  SLAVE2 


C 


CALL    INVERT  (A, F ,R?OSr , RPO STL,RSLAV2, ESLAL2, F AZ2, BAZ 2, 
1    DIST2) 
UDIST2    =     (RN    *    DIST2)     /    UNCOV 
CALL    SECFAC  (UDIST2,SF2) 
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C  COMPOTE    THE    RATES    AND    COMPARE    TO    THE    OBSERVED    RATES 

C 

TDM1    =    UDIST1    -    ODISTM    +    SF1    -    SFM    ■•■    DELAYl 

DIPF1    =    TDM1     -    EAT  El 
C 

TDM2    =    UDIST2    -    UDIST M    +    SF2    -    SFM    +    DSLAY2 

DIFF2    =    TDM2    -    RATE2 
C 

C  COMPUTE    THE    LANE    WIDTH    IN    METERS    BASED    ON    EQUATION 

C  4.20    IN    ELECTRONIC    SURVEYING    AND    NAVIGATION    - 

C  LAURILA,    PAGE    94. 

C 

BR1    =    DABS  (FAZ1    -    FAZ  M) 

3R2    =    DABS(FAZ2    -    FAZ  M) 

WIDTH1    =     (DIFF1     *    UNCOV    *    0.5)     /    DSIN  (BR1    *    0.5) 

WIDTH2    =     {DIFF2    *    UNCOV   *    0.5)     /    DSIN  (BR2    *    0.5) 
C 

C  WRITE   THE    POSITION    OF    VESSEL,    COMPUTED    RATE,    OBSERV^ID 

C  RATE,    AND    THE    DIFF    BETWEEN    THEM. 

C 

WRITE  (6^200)     ILAT,IMIN,RSEC,IL0N,ILMIN,RL3EC,RATS1, 

DIFF1,RATE2,DIFF2 
200         F0RMAT(1X,I3,I3,1X,F6.3,I5,I3,1X^F6.3,F12.2,F8.2, 
1    F12.2,F8.2) 

WRITE(6,201) 
C 

C  XMEAN1    AND    XMEAN2    ARE    THE    MEAN    DIFFS    BETWEEN    THE 

C  COMPUTED    RATS    AND    THE    OBERSEHVED.     XSLAV1    AND    XSLAV2 

C  ARE    THE    STORED    DIFFS. 


C 


XMEAN1  =  XMSAN1  +  DIFF1 

XMEAN2  =  XMEAN2  ■•■  DIFF2 

XMEAN3  =  XMEAN3  *  WIDTHi 

XMEAN4  =  XMEAN4  +  WIDTH2 
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XSLAV1  (I)  =  DIFF1 

XSLAV2(I)  =  DIFF2 
550    CONTINUE 
C 

C      COMPUTE  THE  MIAN  AND  STANDARD  DEVIATION 
C 

XMEAN1  =  XMEAN1  /  IREC 

XMEAN2  =  XMEAN2  /  IREC 

XMEAN3  =  XMEAN3  /  IREC 

XMEAN4  =  XMEAN4  /  IREC 
C 

DO  6  00  I  =  1, IREC, 1 

VAR1     =    7AR1    +     ((XSLA7  1(I)     -    XMSAN1)**2) 
600  7AR2    =    VAR2    +     ((XSLA72(I)     -    XMEAN2)**2) 

C 

VAR1    =    VAR1    /    (IREC    -     1.0) 

VAR2    =    VAR2    /     (IREC    -     1.0) 
C 

SD1    =    DSQRT(VARI) 

SD2    =    DSQRT(VAR2) 
C 

WRITE  (6,201) 

WRITE  (6,210)     XMEAN1,SD1,XMEAN3 

WRITE(6,201) 

WRITE(6,211)     XMSAN2,S D2,XMEANa 

WRITE  (6,202) 

210  F0RMAT(1X,'     SLAVE    #1,     MEAN    =    ',?10.3, 

1  •    STANDARD    DEVIATION    =    •  , 

2  F10.3,'    DISTANCE    IN     METERS    =    ',F10.3) 

211  F0RMAT(1X,'     SLAVE    #2,     MEAN    =    ',F10.3, 

1  '    STANDARD    DEVIATION    =    ', 

2  F10.3,'    DISTANCE    IN     METERS    =    •,F10.3) 
STOP 

END 
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C=====:  =  =  ==  =  ==  ==  =  =  =  =  ========  ==  =  =  =  =  =  =  ====  =  ==  =  =  =  =  =  =  =  =  =  =  ==  =  ==  = 

SUBROUTINE    SECFAC  (UTDIST,  SF) 
C 

c 

C  THIS    ROUTINE    WILL    COMPUTE    THE    SEA    SECONDARY    FACTOR 

C  OTDIST    =    TOTAL    DISTANCE 

C  SF    =    SECONDARY    FACTOR 

C 

IMPLICIT    REAL*8     (A-H,  0-Z) 
C 

c 

C  COEFFICIENTS 

C 

AO   =       129.04323 

A1    =    -0.40758 

A2   =      0.00064576813 
C 

BO   =      2.741282 

B1    =    -0.011402 

B2   =    .00032774815 
C 

IF(UTDIST    .GT.     537)     30    TO    10 

SF    =     (    BO    /    UTDIST)     +    B1    +     (    B2    *    UTDIST) 

GO    TO    20 
10  SF    =     (    AO    /    UTDIST)     •»•     A1    +     (    A2    *    UTDIST) 

C 
20  CONTINUE 

RETURN 

END 

C== ======================== ================================: 

SUBROUTINE  INVSR1  (A,  ? INV, GLAT1 , GLON 1 , GLAT2 , GL0N2  ,FAZ , 
1  BAZ,S) 
C 

c 
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Q   4e«:(c    SOLUTION    OF    THE    GEODETIC    INVERSE    PROBLEM    AFTER 

C   ***    T.VINCENTY    MODIFIED    RAINS  FORD'S    METHOD    WITH    HELMERT'S 

C   ***    ELLIPTICAL    TERMS.       EFFECTIVE    IN    ANY    AZIMUTH    AND    AT 

C   ***    ANY    DISTANCE    SHORT    OF    ANTIPODAL    STANDPOINT/FOREPOINT 

C   ***    MUST    NOT    BE   THE    GEOGRAPHIC    POLE 

C 

C   ***    A    IS    THE    SEMI-MAJOR    AXIS    OF    THE    REFERENCE    ELLIPSOID 

C   ***    FINV    IS    THE    FLATTENING    (NOT    RECIPROCAL)     OF    THE 

C   ***    REFERNECE    ELLIPSOID    LATITUDES    AND    LONGITUDES    IN 

C    ***    RADIANS    POSITIVE    NORTH    AND    EAST    FORWARD    AZIMUTHS    AT 

C   ***    BOTH    POINTS    RETURNED     IN    RADIANS    FROM    NORTH    GEODESIC 

C   ***    DISTANCE    S    RETURNED    IN    UNITS    OF    SEMI-MAJOR    AXIS    A 

C 

C   ***    PROGRAMMED    FOR    CDC-66 00    BY    LCDR    L.PFSIFER    NGS 

C   ***    ROCKVILLE    MD    18FEB75.        MODIFIED    FOR    IBM    SYSTEM    360 

C   ***    BY    JOHN    G    3ERGEN    NGS     ROCKVILLE    MD    7507. 

C 

IMPLICIT    REAL*8     (A-H,0-Z) 

DATA    EPS/0. 5D-1 3/, PI/ 3. 14 1 59265 35  89 8D0/ 

TW0PI=2.*PI 

R=1.-FINV 

TU1  =  R*DSIN(GLAT1)  /DCOS  (GLATI) 

TU2=R*DSIN  (GLAT2) /DCOS  (GLAT2) 

CU1=1 ./D3QRT(TU1*TU1+  1.) 

SD1=CU1*Tai 

CU2=1  ./DSQRT(TU2*TU2+ 1.) 

S=CU1*CU2 

3AZ=S*TU2 

FAZ=BAZ*TU1 

X=GL0N2-GL0N1 
100    SX=DSIN  (X) 

CX=DC0S  (X) 

TU1=CU2*SX 

TU2=BAZ-SU1*CU2*CX 
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SY=DSQRT(Tai*T0  1+TU2*TU2) 
CY=S*CX+FAZ 
Y  =  DATAN2(SY,CY) 
SA=S*SX/SY 
C2A=-SA*SA-H. 
CZ=FAZ+FAZ 

IF(C2A.GT.O.)  CZ=-CZ/C2A+CY 
E=CZ*CZ*2.-1. 

C=  ((-3.*C2A  +  4  .)  *FINV4-a.)  *C2A*FINV/16 . 
D  =  X 

X=  ((E*CY*C  +  CZ)  *SY*C-«-Y)  *SA 
X=  (1  .-C)  *X*FINV+GL0N2-GL0N1 
IF  (DABS  (D-X)  .GT.EPS)     GOTO    100 
FAZ=DATAN2(T01,TU2) 
IF  (FAZ.GE.TWOPI)     FAZ= FAZ-T WOPI 
IF  (FAZ.LT.O.DO)     FA2=F AZ+TWOPI 
BAZ=DATAN2(CU1*SX,BAZ*CX-S01*Ca2)  +PI 
IF  (BAZ.GE.THOPI)     B AZ= BAZ-T WOPI 
IF  (BAZ.LT.O.DO)     BAZ=B  AZ+IWOPI 
X  =  DSQRT  (  (1./R/R-1.)  *C2A+1.)  +1. 
X=  (X-2.)/X 
C=1.-X 

C=  (X*X/U.  +  1.) /c 
D=(0.375*X*X-1.  )  *X 
X=E*CY 
S=1.-E-E 

S=  ((  (  (SY*SY*4.-3.)  *S*CZ*D/6.-X)  *D/4.+CZ)  *SY*D  +  Y)  *C*A*R 
RETURN 
END 
SENTRY 
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DATA    SET    EXAMPLE    -    PEOGRAM    LORAN 

SHIP'S  POSITION               OBSERVED    LORAN    RATES    &    ERRORS     (E) 

LATITUDE  LONGITUDE              9940-Y            EY  9940-W            EW 

(D-M-S)  (D-M-S)                 ()1S6C)  (;isec)        (;isec)  (usee) 

36    U3    45.800  121    55    27.160       42789.31  -0.49  16294.04  -1.06 

36    44      3.400  121    55    32.340       42791.13  -0.38  16293.46  -1.10 

36    44    21.180  121    55    37.390       42793.04  -0.38  16292.73  -0.99 

36    44    37.490  121    55    46.950       42795.13  -0.58  16292.03  -1.06 

36    44    53.260  121    55    57.710  42796.93  -0.51  16291.43  -1.27 
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APPENDIX    C 
PROGRAM    LOPLC 

C  PROGRAM    LOPLC 

C 

C  PROGRAM    COMPUTES    LINE    OF    POSITIONS    FOR    LORAN-C    USING 

C  THE    SECONDARY    FACTOR      (SF)     AND    ADDITIONAL    SECONDARY 

C  FACTOR     (ASF)     THE    SF    AND    ASF    ARE    BASED    ON 

C  SEMI-EMPIRICAL    TD    GRID. 

C 

IMPLICIT    REAL*8     (A-H,0-Z) 

DIMENSION    PHI  (100)  ,ELON  (100)  ,DISTM  (100)  ,DIST1  (100)  , 
1    DIST2(10  0)  ,AXIS  (13)  ,  RF  (13)  ,AZ1  (100)  ,  AZ2  (1  00)  ,  AZM  (1  00) 

DATA    RHOSEC,PI, UNCOV/2. 06 264806  2U7 1D05, 
1    3.1U15926535898D0,299.693D0/ 
C 
cccccc* ******  ********  ******  ************* 

DATA    AXIS/6. 3782064D06, 6. 378388         D06, 6. 377 397 1 55D06 , 

1  6.37816    D06,6. 37816D06^6.3782U91U5D06,6.378165D06, 

2  6.378166D06,6. 378165    D06 , 6. 378 1 a5D06 , 6 . 377563UD06, 

3  6.378245D06,6. 3781350D06/ 

DATA    RF/6. 3565838  D06,2.97         D 02, 2 . 99 15 28 1285D02 , 

1  2.98  25D0  2,2. 98 247 167 4  27D 02, 2 .93 465D02 , 2 .9825D02 , 

2  2.983    D02,2.983D02,2 .9825    D02, 6 , 3562569D06 , 

3  2.983    D02,2.98  26D02/ 
C*****ELLIPSOID    OPTION    NUMBER 

8.  MERCURY 

9.  MARSHALL  ISLAND 

10.  NAVY  8D 

11.  AIRY 

12.  KRASSOWSKI  19U0 

13.  WGS  1972 

88 


C 

1. 

CLARKE    1866 

c 

2. 

INTERNATIONAL    (HAYFORD) 

c 

3. 

BESSEL    1841 

c 

4. 

AND     (AUSTRALIAN) 

c 

5. 

1967    REFERENCE 

c 

6. 

CLARKE    1830    MOD 

C  7.    SAO 

C 

cccccc* ************ ********  ************* 

C  CC1-2    =    ELLIPSOID    NUMBER     (K) 

C  CC3-5    =    NUMBER    OF    POINTS    ALONG    COAST     (IREC) 

C 

READ(5,100)     K,IREC,ISTA1,ISTA2 

100  FOEMAT(1X,I2,I4  ,213) 
TW0PI=2.*PI 

A  =  AXIS(K) 

F=1./RF  (K) 

IF  (F.LT-3.D-3)     F=(A-1./F)/A 
C 

C  READ    MASTER    STATION    POSITIONS. 

C 

READ  (5,  105)     ILATM,IMINM,RSECM,ILONM,ILMINM,RLSECM 
105         FORMAT (1X,I3, 1X,I2,1X,F5.2,1X,I4,1X,I2,1X,F5.2) 

RMAST    =     ( (IA3S (ILATM)     *    60    +    IMINM)     *    60    +    RSSCM)     / 
1    RHOSEC 

IF  (ILATM    .LT.    0)     RMAST    =    -RMAST 

RMASTL    =     ((TABS  (ILONM)     *    60    +    ILMINM)     *60    +RLSECM)     / 
1    RHOSEC 

IF(ILONM    .GT.    0)     RMASTL   =    TWOPI    -    RMASTL 
C 

C  READ    THE    FIRST    SLAVE    STATION    AND    CHANGS    THE    LATITUDE 

C  AND    LONGITUDE    INTO    RA DI ANS . DELAY    IS    THE    CODING    DELAY 

C  AND    BIAS    IS    THE    OFFSET    IN    MICROSECONDS. 

C 

READ  (5, 101)     ILATS,IMINS,RSSCS,ILONS,ILMINS,RLSECS, 
1    DELAY1,3IAS1 

101  FCRMAT(1X,I3, 1X,I2,1X,F5.2, 1X,I4,  1X,I2,1X,F5.2, 
1    ?9.2,F7.3) 

RSLA71    =     ( (lABS  (ILATS)     *    60    +    IMINS)     *    60    *    RSSCS)     / 
1    RHOSEC 
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IF(ILATS    .LT.    0)     RSLAV1    =    -RSLAV1 

RSLAL1    =     (  (lABS  (ILONS)     *    60    +    ILMINS)     *    60    +    RLSECS)     / 
1    RHOSEC 

IF  (ILONS    .GT.    0)    RSLAL1    =    TWOPI    -    RSLAL1 
C 

C  READ    THE    SECOND    SLAVE    STATION    AND    CHANGE   THE    LATITUDE 

C  AND    LONGITUDE    INTO    RADIANS. 

C 

READ  (5, 101)     ILATS,IMINS,RSECS, ILONS, ILMINS, RLSECS, 
1    DELAY2,BIAS2 

RSLAV2    =     ((TABS  (ILATS)     *    60    +    IMINS)     *    60    +    RSECS)     / 
1    RHOSEC 

IF(ILATS    .LT.    0)     RSLAV2    =    -RSLAV2 

RSLAL2    =     ( (lABS (ILONS)     *    60    +    ILMINS)     *    60     +    RLSECS)     / 
1    RHOSEC 

IF  (ILONS    .GT.    0)     RSLAL2    =    TWOPI    -    RSLAL2 
C 

C  READ    COAST    POINT   LAT     AND    LONG    AND    CONVERT    TO    RADIANS. 

C  STORE   LAT    IN    PHI    AND     LONG    IN    ELON. 

C 

DO    5  00    I=1,IRZC,1 

READ     (5,102)     INUM,ILAT,IMIN,RSSC,ILON,ILMIN, RLSSC 
102  FORM  AT ( 1 X, 14, 14, 13, F7. 3, 1  5, 13, F 7. 3) 

PHI  (I)     =     ((lABS(ILAr)      *    60    -<■    I2IIN)     *    60    +    RSSC)     / 
1    RHOSEC 

IF(ILAT    .LT.     0)     PHI  (I)     =    -PHI  (I) 

ELON  (I)     =     ((IABS(IL3N)     *    60    +    ILMIN)     *    60    ^-    RLSEC)     / 
1    RHOSEC 

IF(ILON    .GT.     0)     ELON  (I)     =    TWOPI    -    ELON  (I) 
C 

C  COMPUTE    DISTANCE    FROM    COAST    POINT    TO    MASTER    AND 

C  SLAVE    STATIONS.     SET    UP    COMPUTATION    AND    CALL    INVER1 


C 


?1    =    ?HI(I) 

90 


El  =  ELON(I) 

CALL  INVER1  (A, F, Pi, SI  ,RMAST,R»1ASTL  ,  AZF  ,  AZB  ,  S) 

DISTM(I)  =  S 

AZM(I)  =  AZF 

CALL  INVER1 (A,F,P1,E1 ,RSLAV 1 , RSLAL1 , AZF, AZB , S) 

DIST1  (I)  =  S 

AZ1(I)  =  AZF 

CALL  INVER1 (A,F,P1,E1 ,RSLA72, RSLAL2 , AZF , AZB , S) 

DIST2  (I)  =  S 

AZ2(I)  =  AZF 

WRITE  (6,210)  I,INUM,ILAT,IMIN,RSEC,IICN,ILi5IN,  RLSEC, 
1  AZM (I)  ,AZ1  (I)  , AZ2  (I) 
219    FORM  AT (IX, 14, 14, 14, 1 3 ,F7. 3 , 15 , 13 , F7. 3 , 3F10 . 3) 
500    CONTINUE 
C 

C      READ  THE  NUMBER  OF  RECORDS  AND  THEIR  POSITIONS. 
C      COMPUTE  THE  TOTAL  DISTANCE  OVER  LAND  AND  SEA. 
C       TDISTM  =  TOTAL  DISTANCE  FROJl  FASTER  TO  DATA  POINT. 
C       TDISTS  =  TOTAL  DISTANCE  FROM  SLAVE  TO  DATA  POINT. 
C 

C      COMPUTE  THE  FORWARD  AND  BACK  AZIMUTHS 

C        FAZM  AND  FAZS  =  FORWARD  AZIMUTH  TO  THE  MASTER  AND 
C        SLAVE  STATION. 

C        EAZM  AND  3AZS  =  BACK  AZIMUTH  TO  THE  MASTER  AND 
C        SLAVE  STATION. 
C 

READ  (5,  103)  JREC 

103  FORM  AT (IX, 14) 
C 

DO  550  1=1, JREC, 1 

READ  (5,  104)  ILAT  ,IilIN,  RS  EC  ,  HON,  ILMIN  ,  RLSEC  ,  RATS  1 ,  HAT  E2 

104  FORM  AT (IX, 13,13,  F6. 3,  15,13, F6. 3, F 9.3, F9.  3) 
C 

C  CONVERT    LAT    AND    LONG     TO    RADIANS. 
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RPOST    =     ((lABS(ILAT)      *    60    +    IMIN)     *    60    +    RSEC)     / 
1    RHOSEC 

IF(ILAT    .LT.     0)     RPOSr    =    -RPOST 

BPOSTL    =     ((lABS(ILON)     *    60    +    ILMIN)     *    60    +    RLSEC)     / 
1    RHOSEC 

IF(ILON    .GT.     0)     RPOSTL    =    TWOPI    -    RPOSTL 
C 

C  COMPUTE    DISTANCES    AND    AZIMUTHS    FROM    THE    OBSERVED 

C  POINT    TO    THE    INTERPOLATED    SHORE    POINT    AND    WRITE. 

C 

ISTA    =    1 

CALL    SHORPT (I  RE C^RPOS T, RPOSTL, R MAST, RMASTL , UTDISM, 
1    USDIST,ULDIST,FAZM,B  AZM) 

CALL    SECFAC(UTDISM,aSDIST,ULDI3T,BAZM,ISTA,SFM) 
C 

CALL    SHORPT (I REC,RPOS T, RP0STL,RSLAV1 , RSLAL1 ,UTDIS1, 
1    USDIST,ULDIST,FAZ1,B  AZ1) 

CALL    SECFAC (UTDIS 1 ,as DIST , ULDIST, BAZ1 ,ISTA 1 ,SF1) 
C 

CALL    SHORPT (I  RE C,RPOS T, RPOSTL, RSL A72 , RSLAL2 , UTDIS2, 
1    USDIST,ULDIST,FAZ2,B  AZ2) 

CALL    SECFAC (UTDIS2 ,as DIST , ULDIST, BAZ2 ,ISTA2 , SF2) 
C 

C  COMPUTE    THE    RATES    AND    COMPARE    TO    THE    OBSERVED    RATES 

C 

TDM1    =    UTDISI    -    UTDISM    +    SF1    -    SFM    +    DELAY1    +    BIASI 

DIFF1    =    TDMl    -    RATE1 
C 

TDM2    =    UTDIS2    -    UTDISM    +    SF2    -    SFM    +    DELAY2    +    EIAS2 

DIFF2    =    TDM2    -    RATE2 
C 

C  COMPUTE    THE    LANE    WIDTH    IN    METERS    BASED    ON    EQUATION 

C  4.20    IN    ELECTRONIC    SURVEYING    AND    NAVIGATION     -    LAUEILA, 
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C  PAGE    94. 

C 

BR1    =    DABS  (FAZ1    -    FAZ  M) 

IF(BR1     ,GT.     PI)     BRI    =    TWOPI    -    BR1 

BR2    =    DABS  (FAZ2    -    FAZ  M) 

IF(BR2    .GT.     PI)     BR2    =    TWOPI    -    BR2 

WIDTH1    =     (DIFF1    *    UNCOV   *    0.5)     /    DSIN(BR1    *    0.5) 

WIDTH2    =     (DIFF2    *    UNCOV    *    0.5)     /    DSIN  (BR2    *    0.5) 
C 

C  WRITE   THE    POSITION    OF    VESSEL,       OBSERVED    RATES, 

C  AND    THE    DIFFERENCES    BETWEEN    THEM. 

C 

WRITE  (6,200)     ILAT,IMIN,RSEC,IL0N,ILMIN,P.LSEC,RATS1, 
1    DIFF1,RATE2, DIFF2 
200         F0RMAT(1X,I2,I3,F7.3,  Iit,l3,F7.3,?10.2,F6.2,F10.2,F6.2) 
C 

C  XMEAN1    AND    XMEAN2    ARE    THE    MEAN    DIFFS    BETWEEN    THE 

C  COMPUTED    RATE    AND    THE    OBERSERVED.     XSLAV1    AND    XSLAV2 

C  ARE    THE    STORED    DIFFS, 

C 

XM2AN1    =    XMEAN1    +    DIFF1 

XMEAN2    =    XMEAN2    +    DIFF2 

XMEAN3    =    XMEAN3    +    WIDTH  1 

XMEANU    =    XMEANU    +    WIDTH2 

XSLAV1  (I)     =    DIFF1 

XSLAV2(I)     =    DIFF2 
550         CONTINUE 
C 
C  COMPUTE    THE    MEAN    AND     STANDARD    DEVIATION 


C 


XMEAN1  =  XMEAN1  /  JREC 

XMSAN2  =  XMEAN2  /  JREC 

XMEAN3  =  XMEAN3  /  JREC 

XMEAN4  =  XMEANU  /  JREC 
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c 

DO    600    I    =    1, JREC, 1 

VAR1  =  VAR1  +  ((XSLA7  1(I)  -  XMEAN1)**2) 
600  VAR2  =  VAR2  +  ((XSLAV2(I)  -  XMEAN2)**2) 
C 

7AR1    =    VAR1    /    (JREC    -     1.0) 

VAR2    =    VAR2    /    {JREC   -     1.0) 
C 

SD1    =    DSQRT  (VAR  1) 

SD2    =    DSQRT  (VAR2) 
C 

WRITE(6,201) 

WRITE  (6,210)     XMEAN1,SD1,XMEAN3 

WRITE  (6,201) 

WRITE  (6,211)     XMEAN2,SD2,XMEANa 

WRITE  (6,202) 

210  F0RMAT(1X,'     SLAVE    #1,     MEAN    =    •,F10.3, 

1  •    STANDARD    DEVIATION    =    •, 

2  F10.3,'     DISTANCE    IN    METERS    =    »,F10.3) 

211  F0RMAT(1X,'     SLAVE    #2,     MEAN    =    «,F10.3, 

1  '    STANDARD    DEVIATION    =    •, 

2  F10.3,'    DISTANCE    IN     METERS    =    •,?10.3) 
WRITE(6,201) 

201  FORMAT  (ia    ) 

202  FORMAT(IHI) 
STOP 

END 
C==== ======================  =======  ===  =  ====  =====  =  ====  =  ======== 

SUBROUTINE    SHORPT (IRE C, RPOST , RPOSTL, RCOUT, RCONTL , UDIST , 

1    USD,ULD,FA,3A) 
Q 

C 

C  SUBROUTINE    SHORPT    WILL    SELECT    A    POINT    ALONG    THE    SHORE 

C  WHICH    IS    OUTLINED    FROM    NORTH    BY    SELECTED    POINTS    WITH 
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C  KNOWN    LATITUDES    AND    LONGITUDES.        THE    SHORE    POINT    IS 

C  INTERPOLATED    BETWEEN    TWO    KNOWN    POINTS    USING    THE    TOTAL 

C  DISTANCE    BETWEEN    THE    POSITION    AND    THE    CONTROL    STATION 

C  AND    THE    AZIMUTH    BETWEEN    THE    SHORE    POINTS    AND    THE 

C  RECEIVERS    POSITION. 

C 

IMPLICIT    REAL*8     (A-H,0-Z) 

DIMENSION    FAZM(IOO) 

COMMON/SHORE/PHI  (100)   ,  ELON  (100)  ,  UNCOV,RN,  A  ,  F 

DATA    PI/3,  1U159  26535398D0/ 
C 

TWOPI    =    2.0    *    PI 
C 

CALL    INVER1  (A,F ,RPOST ,RP0STL,R30NT, RCONTL, FA , B A, RDISTT) 

UDIST    =     (RN    *    RDISTT)     /   UNCOV 
C 

RCOMP    =    99999.99 
C 

DO    10    J=1,IREC, 1 

P1    =    PHI (J) 

El    =    ELON  (J) 

CALL    INVER1  (A, F, PI, El  ,RCONT, 3C3NTL , AZF, A23, RDISTL) 

CALL    INVER1  (A,F  ,RPOST ,RP0STL,P1 ,E1 , FAZ, BAZ , RDISTS) 

USD    =     (RN    *    RDISTS)     /    UNCOV 

ULD    =     (RN    *    RDISTL)     /    UNCOV 

FAZM  (J)     =    AZF 
C 

C      COMPUTE  THE  DIFFERENCE  BETWEEN  THE  TOTAL  DISTANCE 
C      (UDIST)  AND  THE  SUMMATION  OF  THE  DISTANCE  OVER  THE  WATER 
C      (USD)  AND  THE  DISTANCE  OVER  THE  LAND  (ULD)  .   IT  THE 
C      DISTANCE  IS  LESS  THAN  RCOMP,  UPDAT  RCOMP  AND  JSTA. 
C      JSTA  IS  THE  CLOSEST  POINT  ALONG  THE  SHORELINE  WHICH  IS 
C      NEAR  THE  EM  PROPAGATION  PATH. 
C 
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FDIFP    =    DABS(UDIST    -      (USD    +    ULD)  ) 

IF(FDIFF    .GT.    RCOMP)      GO    TO    10 

RCOMP    =    FDIFF 

JSTA  =  J 
10  CONTINUE 
C 

C      NOW  DETERMINE  THE  INTERPOLATED  LATITUDE  AND  LONGITUDE 
C      SHORE  POINT  USING  AZIMUTH  PERCENTAGE. 


C 


TOPPER  =  JSTA  -  1 

ILOWER  =  JSTA  ■«■  1 

AZMU    =    FAZM(IUPPER) 

AZML    =    FA ZM  (ILOWER) 

IF  (AZMU    .GT.     FA    -AND.     FA    .GE.    FAZM(JSTA))     ICH    =    lUPPER 

IF  (FAZM  (JSTA)     .  GE.    FA     .AND.    FA    .GT.     AZML)     ICH    =    ILOWER 

RADJ    =    1.00    -    DABS  ((FAZM(ICH)     -    FA)     /     (FAZM  (ICH)     - 
1    FAZM  (JSTA)  )  ) 

RNWPHI    =    PHI  (JSTA)     ^     (  (PHI  (ICH)     -    PHI  (JSTA)  )     *    RADJ) 

IF  (ELON  (JSTA)     .  LE.    ELON(ICH))     RADJ    =    1.0    -    RADJ 

RNWELN    =    ELON  (JSTA)     +     ((ELON  (ICH)     -    ELON  (JSTA))     *    RADJ) 

CALL    TODMS (RNWPHI, IDG  ,MIN, SEC) 

RHOLD    =    TWOPI    -    RNWELN 

CALL    TODMS (RHOLD, IDGL,MINL,SSCL) 

CALL    INVER1  (A,F  ,RNWPH I, RNWELN, RCONT, RCONTL , AZ, BZ , RDISTL ) 

CALL    INVER1 (A,F,RPOST , RPOSTL,RN WPHI, RNWELN , AZ, 3Z, RDISTS ) 

USD    =     (RN    *    RDISTS)     /    UNCOV 

ULD    =     (RN    *    RDISTL)     /    UNCOV 

RETURN 

END 
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SUBROUTINE    SECF AC(UTD 1ST,  USDIST , ULDIST, AZI , I STA , SF) 
C 

c 

C  THIS  ROUTINE  WILL  COMPUTE  THE  LAND/SEA  SECONDARY  FACTOR 

C  UTDIST  =  TOTAL  DISTANCE 

C  USDIST  =  DISTANCE  OVER  THE  SEA  WATER  PATH 

C  ULDIST  =  DISTANCE  OVER  THE  LAND  PATH 

C  AZI  =  AZIMUTH  FROM  NORTH. 

C 

C  MASTER  =  1  (ISTA) 

C  W       =2 

C  X       =3 

C  Y       =4 

C 

IMPLICIT  REAL*8  (A-H,0-Z) 

51  =     (.795    /    USDIST)      +    0.439    +     (.00245    *    USDIST) 
C 

IF  (UTDIST    .GT.     540)     GO    TO     10 

52  =     (3.188    /    UTDIST)     -    0.594    +     (.000329    *    UTDIST) 
GO    TO    20 

10  52   =     (128.8    /   UTDIST)     +    0.187    +     (.000652    -    UTDIST) 

C 

20  IF(USDIST    .GT.     540)     GO   TO    30 

53  =     (3.188    /   USDIST)     -    0.594    +     (.000329    *    USDIST) 
GO   TO    40 

30  S3    =     (128.8    /   USDIST)     +    0.187    +     (.000652    *    USDIST) 

C 

40  IF  (ULDIST    .GT.     540)     GO    TO    50 

54  =     (3.188    /   ULDIST)     -    0.594    +     (.000329    *    ULDIST) 
GO    TO    60 

50  S4    =  (128.8    /   ULDIST)  +    0.187    *■     (.000652    *    ULDIST) 

C 

60  S5   =  1.428    +     (.00158  *    UTDIST) 

S6    =  1.428    -••     (.00158  *    ULDIST) 
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TAZI    =    2-0    *    AZI 

GOTO    (70,80,90,100) , ISTA 
70  RHOLD    =  (1.010*DSIN  (AZI) )    -     ( . 1 96*DCOS  (AZI) ) 

1  -     (.893*DSIN(TAZI) )     -     ( . 355*DCOS (TAZI) ) 

GO   TO    200 
80  RHOLD    =  (.323*DCOS  (AZI)  )     -     (.  711  *DSIN  (TAZI)  ) 

GO   TO    200 
90  RHOLD    =  (.9a2*DCOS (TAZI) ) 

GO    TO    200 
100         RHOLD    =  (.588*DSIN  (TAZI)  ) 

C 
200         S5    =    S5    +    RHOLD 

S6    =    S6    +    RHOLD 

SF    =    0.5    *     (S5    ^    S6    -    SI    +    S2    +    S3    -    SU) 

WRITZ  (6,5  00)     S1 ,S2,S3  ,SU,S4,S5,S6 ,S? 
500         FORMAT(1X,7F15.  5) 

RETURN 

END 
C== ============================================== ==========^ 

SUBROUTINE  INVERT  ( A, ? INV, GLAT1 , GLON 1 ,GLAT2 , GL0N2  ,FAZ, 

1  3AZ,S) 
C 

See   Appendix   B    for   subroutine   INVSR1. 
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short   points   around   Monterey  Bay,    California.      The   point 
are    used  to   interpolate   geodetic   point    for   computation 

of   the  distance   over   land   and  the   distance   ov^r   sea. 

NO.         LATITUDE  LONGITUDE 

1  36    57    18.606  122    05  37.525 

2  36    56    59.264  122    03  01.817 

3  36    57    05.076  122    01  31.701 
U  36    57    49.538  122    01  07.857 

5  36    57    17.949  121    58  19.830 

6  36    58    08.589  121    57  07.288 

7  36    58    32.140  121    55  10. 083 

8  36    58    01.498  121    53  57.390 

9  36    56    46.115  121    52  22.313 

10  36    55    38.140  121    51  24.399 

11  36    53    13.806  121    49  46.743 

12  36    49    38.384  121    47  48.395 

13  36    47    39.241  121    47  10.818 

14  36    46    27.554  121    47  39.637 

15  36    44    56.717  121    47  52.416 

16  36    41     14.439  121    48  32.642" 

17  36    39    17.211  121    49  28.533 

18  36    37    31.128  121    50  31.723 

19  36    36    23.446  121    51  34.833 

20  36    36    03.628  121    52  50.879 

21  36    36    24.782  121    53  48.453 

22  36    37    18.151  121    54  11.628 

23  36    38    00.300  121    55  57.538 
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DATA    SET    EXAMPLE    -    PROGRAM  LOPLC 

SHIP'S    POSITION  OBSERVED    LORAN    RATES    &    ERRORS     (E) 

LATITUDE              LONGITUDE  9940-Y  EY  9940-W            EM 

(D-M-S)                    (D-M-S)  (;isec)  (;isec)  (jjusec)  (jasec) 

36    43    45-800    121    55    27.160  42789.34  1.14  16294.04  1.32 

36    44      3.400    121    55    32.340  42791.13  1.25  16293.46  1.29 

36    44    21.180    121    55    37.390  42793.04  1.25  16292.73  1.40 

36    44    37.490    121    55    46.950  42795.13  1.05  16292.03  1.33 

36    44    53.260    121    55    57.710  42796.93  1.12  16291.43  1.11 
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APPENDIX    D 
PROGRAM   LORTAB 

C  PROGRAM    LORTAB 

C 

C  PROGRAM    COMPUTES    LINE    OF    POSITION    FOR    LORAN-C    USING 

C  SF    SALT    WATER    CORRECTION    FACTOR     AND    DMAHTC    CALCULATED 

C  OR    FIELD    OBSERVED    ASF    CORRECTIONS. 

C 

C 

IMPLICIT    REAL*8     (A-H,  0-Z) 

DIMENSION    AXIS  (13)  ,RF  (13) 

DIMENSION    XSLAV1  (1000)  ,XSLAV2  (1000) 

DATA    RHOSEC,PI, UNCOV,  RN/2 . 06264 8062471  DOS , 
1    3. 14 15 926535 89 800,29  9. 79 24 58D0 ,1 .00 03 3 8 DO/ 

DATA    XMEAN1,XMEAN2,VAR1//AR2/0. OODO, O.OODO, 
1    O.OODO,OO.ODO/ 

DATA    XMEAN3,XMEAN4/0. OODO,O.OODO/ 
C 
CCCCCC*  ************  ********  ♦***********:* 

DATA    AXIS/6. 3782064006,6. 373338  006,6.377397155006, 

1  6.37816    006,6.37816006,6.378249145006,6.373165006, 

2  6,378166006,6.378165    006,6.37814  5006,6.3775634006, 

3  6.378245006,6,3781350006/ 

DATA    RF/6. 3565838  006,2.97         002,2.99152312  35002, 

1  2.9  8  2500  2,2.98  24716  7  4  270  02,2.9  3  46  5002,2.9  8  2500  2, 

2  2.983    002,2.983002,2.9825    002,6.3562569006, 

3  2,983    002,2.9826002/ 
C*****ELLI?SOID    OPTION    NUMBER 

C  1.    CLARKE    1866  8.     MERCURY 

C  2.    INTERNATIONAL    (HAIFORD)  9.     MARSHALL    ISLAND 

C  3.    BESSEL    1841  10.     NAVY    80 
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c 

4.    AND    (AUSTRALIAN)                          11.     AIRY 

c 

5.     1967    REFERENCE                               12.     KRASSOWSKI     1940 

c 

6.    CLARKE    1880    MOD                            13.     WGS    1972 

c 

7.    SAO 

c 

cccccc* ************ ********  ************* 

1        c 

CC1-2    =    ELLIPSOID    NUMBER    (K) 

c 

CC3-5    =    NUMBER    OF    POINTS    ALONG    COAST     (IREC) 

c 

READ(5,100)     K,IREC,ID1,ID2 
100  F0RMAT(1X,I2,I4,2I3) 

TW0PI=2.*PI 

A=AXIS(K) 

F=1./RF  (K) 

IF  (F.LT.3.D-3)     F=(A-1./F)/A 
C 

C  READ    MASTER    AND    SLAVE    STATIONS    POSITIONS 

C  THE    FIRST    RECORD    IS    THE    NUMBER    OF    MASTER    AND    SLAVS 

C  STATIONS    FOR    THE   PARTICULAR    CHAIN. 

C 

WRITS  (6,202) 
202         FORMAT(IHI) 

WRITE(6,201). 
201  F0RMAT(1H    ) 

READ  (5, 105)     ILATM,IMINM,RSSCM,ILONM,ILMINM,RLSECM 

WRITE  (6,105)     ILATM,IMINM,RSSCM,ILONM,ILMINM,HLSECM 
105         F0RMAT(1X,I3, 1X,I2,1X,F5. 2, 1X,I4, 1X,I2,1X,F5.2) 

RMAST    =     (  (lABS(ILATM)      *   60    +    IMINM)     *    60    +    RSECM)     / 
1    RHOSEC 

IF(ILATM    .LT.    0)     RMAST    =    -RMAST 

RMASTL    =     ((TABS  (ILONM)     *    60    +    ILMINM)     *60    +RLSSCM)     / 
1    RHOSEC 

IF  (ILONM    .GT.    0)     RMASTL   =    TWOPI    -    RMASTL 
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C      READ  THE  FIRST  SLAVE  STATION  AND  CHANGE  THE  LATITUDE 

C      AND  LONGITUDE  INTO  RADIANS. 

C 

READ  (5,  101)     IIATS,IMINS,ESECS,ILONS,ILMINS,RLSECS, 
1    DELAY1 

WRITE  (6,1 01)     IIATS,IMINS,RSECS,ILONS,ILMINS,RLSECS, 
1    DELAY1 
101         FORMAT (IX, 13, IX ,12 , 1 X ,F5. 2, IX, I  4 , 1 X,I2, IX, F5 . 2 , F9 . 2) 
RSLA71    =     (  (lABS  (ILATS)     *    60    +    IHINS)     *    60    +    5SECS)     / 
1    RHOSEC 
IF(ILATS    .LT.    0)    RSLAV1    =    -PSLAV1 

RSLAL1    =     ((lABS  (ILONS)     *    60    +    ILMINS)     *    60    +    RLSECS)     / 
1    RHOSEC 
IF  (ILONS    .GT.    0)     RSLAL1    =    TWOPI    -    R3LAL1 
C 

C  READ    THE    SECOND    SLAVE    STATION    AND    CHANGE    THE    LATITUDE 

C  AND      LONGITUDE    INTO    RADIANS. 

C     - 

READ  (5, 101)     ILATS,IMINS,RSECS, ILONS, ILMINS, RLSECS, 
1    DELAY2 

WRITE  (6,101)     ILATS, iaiNS,RSECS, ILONS, ILMINS, RLSECS, 
1    DELAY2 

RSLAV2    =     ((lABS  (ILATS)     *    60    +    IMINS)     *    60    +    RSECS)     / 
1    RHOSEC 
IF  (ILATS    .LT.    0)     RSLAV2    =    -RSLAV2 

RSLAL2    =     ((IA3S  (ILONS)     *    60    +    ILMINS)     *    60    +    RLSECS)     / 
1    RHOSEC 
IF  (ILONS    .GT.    0)     RSLAL2    =    TWOPI    -    RSLAL2 
WRITE  (6,201) 
C 

C  READ    THE    RECORDS    POSITIONS    AND    OBSERVED    RATES,     THEN 

C  COMPUTE    THE    TOTAL    DISTANCE. 


C 


DO    5  50    I=1,IREC,1 
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READ  (5, 10U)     ILATrlMIN ,RSEC, ILOM , ILMIN, RLSEC , 
1    RATE1,RATE2 
10U         FORM  AT ( IX, 13, 13, F6. 3, 15 ,1 3 ,F6  .  3 , 3X, 2F10 . 2) 
C 

C  CONVERT    LAT    AND    LONG    TO    RADIANS. 

C 

RPOST    =     ((lABS(ILAT)     *    60    +    IMIN)     *    60    +    RSEC)     / 
1    RHOSEC 

IF(ILAT    .LT.     0)     RPOST    =    -RPOST 

RPOSTL    =     ((IA3S(IL0N)     *    60    +    ILMIN)     *    60    +    RLSSC)     / 
1    RHOSEC 

IF(ILON    .GT.     0)     RPOSTL    =    TWOPI    -    RPOSTL 
C 

C  COMPUTE    DISTANCES    AND    AZIMUTHS    FROM    THE    OBSERVED    POINT 

C 

C  MASTER 

C 

CALL    INVER1  (A ,F , RPOST , RPOSTL, RM AST, RMASTL, FAZM ,BAZM, 
1    DISTM) 

UDISTM    =     (RN    *    DISTM)     /   UNCOV 

CALL    SECF AC  (UDISTM, SFM) 
C 

C  SLAVE1 

C 

CALL    INVER1 (A,r , RPOST , RPO STL,RSLA V1 , RSLALI , FAZ 1 , BAZ 1 , 
1    DIST1) 

UDIST1    =     (RN    *    DIST1)     /   UNCOV 

CALL    SECFAC  (UDIST1,SF  1) 
C 
C  SLAVE2 


C 


CALL    INVER1 (A, F , RPOST , RPO STL, RS LA V2, RSLAL2 , F AZ2 , 3A22, 
1    DIST2) 
UDIST2    =     (RN    *    DIST2)     /    iJNCOV 
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CALL    SSCFAC{0DIST2,SF2) 
C 

C  DETERMINE    THE    ADDITIONAL    SECONDARY    CORRECTORS    FROM    THE 

C  LORAN-C    CORRECTION    TABLE    FOR    THE    WEST    COAST    CHAIN    99U0 

C 

CALL    TABLE(RPOST,RPOSTL,ID1,ASF1) 

CALL    TABLE(RPOST,RP0STL,ID2,ASF2) 
C 

C  COMPUTE    THE    HATES    AND    COMPARE    TO    THE    OBSERVED    RATES 

C 

TDM1    =    UDIST1    -    UDI3TM    +    SF1    -    SFM    +    ASF1    +    DELAY1 

DIFF1    =    TDM1    -    RATE1 
C 

TDM2    =    UDIST2    -    UDISTM    +    SF2    -    SFM    +    ASF2    -»•    DELAY2 

DIFF2    =    TDM2    -    RATE2 
C 
C 

C  COMPUTE    THE    LANE    WIDTH    IN    METERS    BASED    ON    EQUATION 

C  4.20    IN    ELECTRONIC    SURVEYING    AND    NAVIGATION    - 

C  LAURILA,    PAGE    94. 

C 

BR1    =    DABS  (FAZ1    -    FAZ  M) 

BR2    =    DABS(FA22    -    FAZ  M) 

WIDTH1    =     (DIFF1     *    UNCOV    *    0.5)     /    DSIN(BR1    *    0.5) 

WIDTH2    =     (DIFF2    *    UNCOV    *    0.5)     /    DSIN(BR2    *    0.5) 
C 

C  WRITE   THE    POSITION    OF    VESSEL,    COMPUTED    RATE,    OBSERVED 

C  RATE,    AND    THE    DIFF    BETWEEN    THEM. 

C 

WRITE  (6,200)     ILAT,IMIN,RSEC,IL0N,ILMIN,RLSEC,RATE1, 

DIFF1,RATE2,DIFF2 
200         F0RMAT(1X,I3,I3,1X,F6 . 3,1 5,13, 1 X, F6. 3 ,F 12. 2 , F8 . 2, 
1    F12.2,F8.2) 

WRITE  (6,201) 
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c 
c 

c 
c 
c 


550 

c 

c 
c 


XMEAN1    AND    XMEAN2  ARE    THE    MEAN    DIFFS    BETWEEN    THE 

COMPUTED    RATE    AND  THE    OBEHSERVED.     XSLAVI    AND    XSLAV2 
ARE    THE    STORED    DIFFS. 

XMEAN1    =    XMEAN1    +  DIFF1 

XMEAN2    =    XMEAN2    ••-  DIFF2 

XMEAN3    =    XMEAN3    +  MIDTH1 

XMEAN4    =    XMEANU    +  WIDTH2 
XSLAV1  (I)     =    DIFF1 
XSLAV2(I)     =    DIFF2 
CONTINUE 

COMPUTE    THE    MEAN    AND    STANDARD    DEVIATION 

XMEAN1    =    XMEAN1    /  IREC 

XMEAN2    =-XMEAN2    /  IREC 

XMEAN3    =    XMEAN3    /  IREC 

XMEAN4    =    XMEANa    /  IREC 


600 
C 


DO    6  00    I    =    1,  IREC,  1 

VAR1     =    VAR1    +     ((XSLA7  1(I)      -    XMEAN1)**2) 

VAR2    =    VAR2    ■••     ((XSLA72(I)     -    XMEAN2)**2) 

VAR1    =    VAR1    /    (IREC    -     1.0) 
VAR2    =    VAR2    /    (IREC   -     1.0) 


SD1    =    DSQRT  (VAR  1) 
SD2    =    DSQRT  (VAR2) 

WRITE  (6,201) 

WRITE(6,210)     XMEAN1,SD1,XMEAN3 

WRITE  (6,201) 

WRITS(6,211)     XMEAN2,SD2,XMEAN4 
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WRIT2(6,202) 

210  F0RMAT(1X,»    SLAVE    #1,     MEAN    =    SFIO.S, 

1  •    STANDARD    DEVIATION    =    «, 

2  P10.3,'    DISTANCE    IN     METERS    =    »,F10.3) 

211  F0RMAT(1X,'    SLAVE    #2,     MEAN    =    ',F10.3, 

1  •    STANDARD    DEVIATION    =    », 

2  F10.3,'    DISTANCE    IN    METERS    =    •,F10.3) 
STOP 

C=========  =  ========== ======  ====  ===  ===  =  =  ===  ======  =  =  =  =  =  =  === 

SUBROUTINE    SECF AC(UTD 1ST, SF) 
C 

C 

C  THIS    ROUTINE    WILL    COMPUTE    THE    SEA    SECONDARY    FACTOR 

C  OTDIST    =    TOTAL    DISTANCE 

C  SF    =    SECONDARY    FACTOR 

C 

IMPLICIT    REAL*8     (A-H,0-Z) 
C 
C 

C      COEFFICIENTS 
C 

AO   =       129.0U323 

A1    =    -0.U0758 

A2   =      0.00064576813 
C 

BO    =      2.741282 

B1    =    -0.011402 

B2    =    .000327748 15 
C 

IF(UTDIST    .GT.     537)     GO    TO     10 

SF    =     (    BO    /    UTDIST)     +    B1    +     (    32    *    UTDIST) 

GO    TO    20 
10  SF    =     (    AO    /    UTDIST)     +    A1    +     (    A2    *    UTDIST) 

C 
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20  CONTINUE 

BETURN 
END 

C=="=  =  =  ==  ==  =  ====  =  =  ==:======  =====  =  =  =  ==  =  =  =  =:==:==  =  =  =  =  =  =  =  =  =  =  =  ====  = 

SUBROUTINE    TABLE(RLAT  ,RLON,ID,  ASF) 

Q 

C 

C  SUBROUTINE    TABLE    SELECTS    THE    PROPER    ASF    CORRECTOR    FROM 

C  THE    LORAN-C    CORRECTION    TABLE    PUBLISHED    BY    THE    DEFENSE 

C  MAPPING    AGENCY. 

C  RLAT    =    POSITION    LATITUDE    IN    SECONDS 

C  RLON    =    POSITION    LONGITUDE    IN    SECONDS 

C  ID         =    LORAN-C    CHAIN   IDENTIFIER 

C  W    =     1 

C  X    =    2 

C  Y    =    3 

C  ASF    =    ADDITIONAL    SECONDARY    FACTORS 

C 

C  THE    FOLLOWING    TABLES     OF   ASF    CORRECTORS    ARE    FOR 

C  MONTEREY    BAY,    CALIFORNIA    -    9940    -W,    -X,    -Y, 

C 

IilPLICIT    REAL*8     (A-H,  0-Z) 

DIMENSION    TA3LEW(3,5)   ,TAB  LEX  (3  ,  5)  ,TABLEY  (3  ,  5) 
C 

DATA    TABLSW/    1. 6D0, 1 - 6D0 , 0. ODO , 

1  1. 6D0,1. 5D0,1.4D0, 

2  1. 5D0,1, 4D0,1 .6D0, 

3  1. 4D0,1.  3D0,1 .5D0, 

4  1. 3D0,0.  ODO, O.ODO/ 
C 

DATA  TABLEX/    -0 . 9D0, -  0. 9D0 ,    O.ODO, 

1  -0.9D0,- 0.9D0,-0.3D0, 

2  -1  ,0D0,- 1.0DO,-1.0DO, 

3  -1 .ODO,- 1.2D0,-1 .1D0, 
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a  -1.1D0,     O.ODO,    O.ODO/ 

c 

DATA    TABLEY/    0. 2D0 , 0.  2D0,0 - ODO, 

1  0. 3D0,0. 3D0,0.3D0, 

2  0. 3D0,0. 2D0,0.4D0, 

3  0. 4D0,0.  3D0,0.6D0, 

4  0. 5D0,0.  ODO, O.ODO/ 

DATA    RHOSEC,PI/2,0626  4306  2a71DD5,3.1U159  26  5  358  98D0/ 
C 

TWOPI    =    PI    *    2. 0 
C 

C  CONVEE    RLAT    AND    RLON    TO   SECONDS 

C 

HLAT    =    ELAT    *    RHOSEC 

HLON    =    TWOPI    -    RLON 

HLON    =    HLON    *    RHOSEC 
C 

C  STARTING    LAT    AND    LONG     FOR    SEARCH 

C  LAT   =    37/00/00.0  LONG    =    122/05/00.0 

C 

c 

C  DETERMINE    THE    AS?    CORRECTOR    FOR    THE    LORAN-C    COMBINATION 

C 

C  LATITUDE 


C 


SLAT    =    133200.0 
SLON    =    439500.0 
RMID    =    300.0 
RDIFF    =    150.0 

DO    10    J=1,5, 1 
SLAT    =    SLAT    -    RMID 
ULAT    =    SLAT    +    RDIFF 
VLAT    =    SLAT    -    RDIF? 
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10     IF  (hLAT  .LT.  DLAT  .AND.  hLAT  .3E.  VLAT)  GO  TO  15 

15     CONTINUE 

C 

C      LONGITUDE 

C 

DO  30  1=1,3,1 

SLON  =  SLON  -  RMID 

ULON  =  SLON  +  RDIFF 

VLON  =  SLON  -  RDIFF 
30     IF  (hLON  .LT.  ULON  .iND.  h  LON  .3T.  VLON)  GO  TO  35 
35     CONTINUE 
C 
C      DETERMINE  ASF  CORRECTOR 


C 


ASF  =0.0 

IF  (ID  .EQ.  1)  ASF  =  rABLEW(I,J) 
IF  (ID  .SQ.  2)  ASF  =  TABLEX(I,J) 
IF  (ID  .EQ.  3)  ASF  =  TABLEY(I,J) 

RETURN 
END 


SUBROUTINE  INVER1 ( A, F INV, GLAT1 , GL0N1 ,GLAT2 , GL0N2  ,FAZ, 
1  BAZ,S) 


See  Appendix  B  for  subroutine  IN7ER1 
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DATA    SET    EXAMPLE    -    PROGRAM    LORTAB 

SHIP'S    POSITION  OBSERVED    LORAN    RATES    &    ERRORS     (E) 

LATITUDE  LONGITODE  9940-Y  EY  9940-W  EW 

(D-M-S)  (D-M-S)  (psec)        (fisec)       {fisec)        (usee) 

36  U3    45.800  121    55  27.  160 

36  44      3.400  121    55  32.340 

36  44    21.  130  121    55  37.390 

36  44    37.490  121    55  46.  950 

36  44    53.  260  121    55  57.710 


42789. 

.34 

-0. 

29 

16294, 

.04 

0. 

34 

42791. 

13 

-0. 

18 

16293 

.46 

0. 

30 

42793, 

04 

-0, 

18 

16292, 

.73 

0. 

41 

42795. 

13 

-0. 

38 

16292, 

.03 

0. 

34 

42796. 

93 

-0. 

31 

16291, 

.43 

0. 

13 
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APPENDIX    E 
PROGRAM    ASFSEL 

C  PROGRAM    ASFSEL 

C 

C  PROGRAM    DETERMINES    OBSERVED    ASF    CORRECTORS    BY    SCANNING 

C  DATA    AT    1     DEGREE    LATITODE    AND    LONGITODE    INTERVALS.    THE 

C  ASF    CORRECTORS    ARE    DETERMINED    3Y    SUBTRACTING    THE 

C  CALCULATED    TD    USING    THE   SEAWATER    SECONDARY    FACTOR    FROM 

C  THE    OBSERVED    TD    RATES. 

C 

IMPLICIT    REAL*8     (A-H,0-Z) 

DIMENSION    AXIS  (  13)  ,RF  (13) 

DIMENSION    ASFCR  1  (16,2  6)  ,ASFCR2(16  ,2  6)  , 
1    IN01  (16,26)  ,IN02(16,  26) 

DATA    RHOSEC,PI, UNCOV,  RN/2, 06264 3062471D05, 
1    3. 1 4  15 9265 35 89 8D0, 29 9.79  245800,1 .00  03  38D0/ 

DATA    XMEAN1,XMEAN2,VAR1  //AR2/0.  0  0D0,0.00D0, 
1    0.00D0,00.0D0/ 

DATA    XMEAN3,XMEAN4/0.  00D0,0.00D0/ 
C 
cccccc* ********************  ************* 

DATA    AXIS/6. 3782064DO 6,6.  378383  D06 , 6. 3773971 55D06 , 

1  6.37816    D06, 6. 37816D06, 6. 378249145006,6. 378165D06, 

2  6.373166006,6.378165    D06 , 6 . 378 1 45D06 ,6 . 3775634006 , 

3  6.378245006,6.3781350006/ 

DATA    RF/6. 3565838  006,2.97         D 02 , 2 .99 152312 85D02 , 

1  2. 9  8 2500 2, 2. 98 247 16 7 4270  02,2.9 3 465002,2.9  8  2500 2, 

2  2.983    002,2.983002,2.9825    002,6.3562569006, 

3  2.983    002,2.9326002/ 
C*****ELLI?SOID    OPTION    NUMBER 

C  1.    CLARKE    1866  8.     MERCURY 
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C  2.  INTERNATIONAL    (HAYFORD)  9.  MARSHALL    ISLAND 

C  3.  BESSEL    1841  10.  NAVY    8D 

C  4.  AND     (AUSTRALIAN)  11.  AIRY 

C  5.  1967    REFERENCE  12.  KRASSCWSKI     194  0 

C  6.  CLARKE    1880    MOD  13.  WGS    1972 

C  7.  SAO 

C 

cccccc* ****** ********  ****** ************* 

C  CCl-2    =    ELLIPSOID    NUMBER     (K) 

C  CC3-5    =    NUMBER    OF    POINTS    ALONG    COAST     (IREC) 

C 

READ(5,100)     K,IREC,ID1,ID2 
100  F0RMAT(1X,I2,I4,2I3) 

TW0PI=2.*PI 

A  =  AXIS(K) 

F=1./RF  (K) 

IF  (F.LT.3.D-3)     F=(A-1./F)/A 
C 

C  READ    MASTER    AND    SLAVE    STATIONS    POSITIONS 

C  THE    FIRST    RECORD    IS    THE    NUMBER    OF    MASTER    AND    SLAVE 

C  STATIONS    FOR    THE   PARTICULAR    CHAIN. 

C 

WRITE  (6,202) 
202         F0RMAT(1H1) 

WRITE  (6,201) 
201  FORMAT (1H    ) 

READ  (5, 105)     ILATM,IMINM,RSECM,ILONM, ILMINM, RLSSCM 

WRITE  (6,105)     ILATM,IMINM,RSSCM,ILONM,ILMINM,RLSECM 
105  FCRMAT(1X,I3, 1X,I2,1X ,F5. 2, 1X,I4, 1X,I2, 1X,?5.2) 

RMAST    =     ( (lABS(ILATM)     *   6 0    >    IMINM)     *    60    +    RSECM)     / 
1    RHOSEC 

IF  (ILATM    .LT.    0)     RMAST    =    -RMAST 

RMASTL    =     (  (TABS  (ILONM)     *    60    +    ILMINM)     *60    +RLSECM)     / 
1    RHOSEC 
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IF(ILONM    .GT-    0)     RMASTL   =    TWOPI    -    RMASTL 
C 

C  READ    THE    FIRST    SLAVE    STATION    AND    CHANGE    THE    LATITUDE 

C  AND    LONGITUDE    INTO    RADIANS. 

C 

READ  (5, 101)     ILATS,IMINS,RSECS,ILONS,ILMINS,RLSECS, 
1    DELAY1 

WRITE  (6,101)     ILATS,IMINS,RSECS,ILONS,ILMINS,RLSECS, 
1    DELAY1 
101         FORM  AT (IX, 13, IX, 12, IX, F5.  2, IX,  14, IX, 12, IX, F5. 2, F9. 2) 
RSLAV1    =     ((lABS  (ILATS)     *    60    +    IMINS)     *    60    +    RSECS)     / 
1    RHOSEC 
IF(ILATS    .LT.    0)    RSLAV1    =    -RSLAV1 

RSLAL1    =     ((lABS  (ILONS)     *    60    +    ILMINS)     *    60    +    RLSECS)     / 
1    RHOSEC 
IF  (ILONS    .GT.    0)     RSLAL1    =    TWOPI     -    RSLAL1 
C 

C  READ    THE    SECOND    SLAVS    STATION    AND    CHANGE    THE    LATITUDE 

C  AND      LONGITUDE    INTO    RADIANS. 

C 

READ  (5,  101)     ILATS,  IMINS, RSECS,  ILONS,  ILMINS,  RLSECS, 
1    DELAY2 

WRITE  (6,101)     ILATS, I ^ INS, RSECS, ILONS, ILMINS, RLSECS, 
1    DELAI2 

RSLAV2    =     (  (TABS  (ILATS)     *    60    +    IMINS)     *    60    +    RSECS)     / 
1    RHOSEC 
IF  (ILATS    .LT.    0)     RSLAV2   =    -RSLAV2 

RSLAL2    =     (  (lABS  (ILONS)     *    60    +    ILMINS)     *    60    +    RLSECS)     / 
1    RHOSEC 
IF  (ILONS    .GT.    0)     RSLAL2    =    TWOPI    -    RSLAL2 
WRITE(6,201) 
C 

C  READ    THE    RECORDS    POSITIONS    AND    OBSERVED    RATES,     THEN 

C  COMPUTE    THE    TOTAL    DISTANCE. 
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DO    55  0    I=1,IREC,1 

READ  (5, 104)     ILAT,IMIN ,RSEC,ILON , ILMIN , RLSEC, 
1    RATE1,RATE2 
104         FOEMAT(1X,I3,I3,F6.3,  15,13 ,F6.  3  ,  3X, 2F1 0. 2) 
C 

C  CONVERT    LAT    AND    LONG    TO   RADIANS. 

C 

RPOST    =     ((IA3S(ILAT)      *    60    +    IMIN)     *    60    +    RSEC)     / 
1    RHOSEC 

IF(ILAT    .LT.     0)     RPOST    =    -RPOST 

RPOSTL    =     ((IABS{ILON)      *   6  0    -»•    ILMIN)     *    60    +    RLSEC)     / 
1    RHOSEC 

IF(ILON    .GT.     0)     RPOSTL    =    TMOPI    -    RPOSTL 
C 

C  COMPUTE    DISTANCES    AND    AZIMUTHS    FROM    THE    OBSERVED    POINT 

C 

C  MASTER 

C 

CALL    INVER1  (A ,F  , RPOST  ,  RPOSTL, RM AST ,RMASTL, FAZM , BAZM , 
1    DISTM) 

UDISTM    =     (RN    *    DISTM)     /    UNCOV 

CALL    SECFAC  (UDISTM,SFM) 
C 

C  SLA7E1 

C 

CALL    INVER1  (A,F  ,RPOST ,RPO STL,RSLA71 , RSLAL1  ,  FAZ1  ,  3AZ  1 , 
1    DIST1) 

UDIST1    =     (RN    *    DIST1)     /    UNCOV 

CALL    SECFAC  (UDIST1,SF  1) 
C 
C  SLAVE2 


C 


CALL    INVER1 ( A, F, RPOST ,RPO STL,RS LAV2, RSLAL2, F AZ2, 3AZ2, 
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1  DIST2) 
0DIST2  =  (RN  *  DIST2)  /  UNCOV 
CALL  SECFAC  (UDIST2,SF2) 
C 

C      DETERMINE  WHICH  LAT  AND  LONG  THE  ASF  CORRECTOR 

C      IS  ASSIGNED  TO. 

C 

CALL  ASSIGN (RPOST, RPO STL, JN1, JN2) 
C 

C  COMPUTE    THE    RATES    AND    COMPARE    TO    THE    OBSERVED    RATES 

C 

TDM1    =    DDIST1    -    UDISTM    +    SFI    -    SFM    +    DELAYI 

DIFF1    =    TDM1     -    RATS1 
C 

TDM2    =    UDIST2    -    UDISTM    +    SF2    -    SFK    +    DELAY2 

DIFF2    =    TDM2    -    RATS2 
C 

C  SUM    THE    DIFFERENCES    TO    THE    MATRIX    AND    COUNT    THE    NUMBER 

C  OF    ASF    CORRECTORS    FOR    EACH    BLOCK    TO    LATTER    DETERMINE 

C  THE    MEAN. 

C 

ASFCR1  (JN1,  JN2)     =    ASFCR  1  (  JN  1,  JN  2)     +    DIFFI 

ASFCR2(JN1,JN2)     =    ASFCR2 ( JN1, JN2)     +    DIFF2 

IN01  (JN1,JN2)  =  IN01  ( JN1,  JN2)  +  1 
550  IN02  (JN1,JN2)  =  IN02  ( JN1,  JN2)  +  1 
C 

C  DETERMINE    THE   MEAN    ASF    CORRECTOR    FOR    EACH    LAT    AND    LONG 

C 

DO    650    1=1, 16,1 

WRITE(6,201) 

DO    600   J=1,26,1 

IF(IN01(I,J)     .EQ.    0)      GO    TO    580 

ASFCR1(I,J)     =    ASFCR1(I,J)      /    IN01(I,J) 
580         IF(IN02(I,J)     .EQ.    0)     GO   TO    590 
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aSFCR2(I,J)     =    ASFCR2(I,J)     /   IN02(I,J) 
590         IF(ASFCR1  (I,J)     .EQ.    0.0)     ASFCR1(I,J)     =    9.99 
600         IF(ASFCR2  (I,J)     .EQ.    0.0)     ASFCR2(I,J)     =    9.99 
650         CONTINUE 
C 

C  WRITE    THE    CORRECTORS     IN    MATRIX    FORMAT 

C 
C 

IF(ID1     .EQ.     1)     WRITE  (6,1000) 

IF(ID1    .EQ.    2)     WRITE(6,1001) 

IF(ID1    .EQ.    3)     WRITE(6,  1002) 

WRITE(6,201) 
C 

WRITE  (6,300)  (  (ASFCR1  (I,J)  ,1  =  1  ,  1  6)  ,  J=1 ,  26) 

DO    700    L=1,5, 1 

WRITE(6,201) 
201  FORMAT  (1H    ) 

700         CONTINUE 
C 

IF(ID2    .EQ.     1)     WRITE  (6,1000) 

IF(ID2    .EQ.    2)     WRITE  (6,100  1) 

IF(ID2    .EQ.    3)     WRITE(6,1002) 

WRITE  (6,201) 

WRITE  (6,300)  (  (ASFCR2  (I,J)  ,1  =  1,16)  ,J  =  1,26) 
300         F0RMAT(16F6.2) 

1000  FORMATC    TABLE    FOR    99U0-W    ») 

1001  FORMATC    TABLE    FOR    99  UO-X     ') 

1002  FORMATC    TABLE    FOR    9940-1     •) 
STOP 

END 
C"==  =====  =  ================  ==  =======  =  =  ===  ======  = 

SUBROUTINE    SECF AC( UTD 1ST, SF) 
C 

C 
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C  THIS    ROUTINE    WILL    COMPUTE    THE    SEA    SECONDARY    FACTOR 

C  UTDIST    =    TOTAL    DISTANCE 

C  SF    =    SECONDARY    FACTOR 

C 

IMPLICIT    REAL*8     (A-H,  0-Z) 
C 
C 

C  COEFFICIENTS 

C 

AO  =   129.04323 

A1  =  -0,40758 

A2   =      0.00064576813 
C 

BO    =      2,741282 

B1    =    -0.011402 

32  =  .00032774815 
C 

IF  (UTDIST  .GT.  537)  GO  TO  10 

SF  =  (  30  /  UTDIST)  +  31  +  (  32  *  UTDIST) 

GO  TO  20 
10     SF  =  (  AO  /  UTDIST)  +  A1  +  (  A2  *  UTDIST) 
C 
20     CONTINUE 

RETURN 

END 
C== ===============  =  ==================  =  =  ===  ===  =  ===  =  ===  =  =  = 

SUBROUTINE  ASSIGN (RLA T, RLON, I, J) 
C 

c 

C  SUBROUTINE  ASSIGN  SELECTS  THE  COLUMN  AND  ROW  FOR 

C  THE  LATITUDE  AND  LONGITUDE  OF  THE  RECORD. 
C         RLAT  =  POSITION  LATITUDE  IN  SECONDS 
C         RLON  =  POSITION  LONGITUDE  IN  SECONDS 
C         I     =  COLUMN 
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C  J  =    ROW 

C 

IMPLICIT    REAL*8     (A-H,  0-Z) 

DATA    RHOSEC,PI/2.0  62  6a806  2U71D0  5, 3 . 1 U 1 59265  35898D0/ 
C 

TWOPI   =    PI    *    2. 0 
C  CONVER    RLAT    AND    RLON    TO   SECONDS 

C 

HLAT    =    RLAT    *    PHOSEC 

HLON    =    TWOPI    -    RLON 

HLON    =    HLON    *    RHOSEC 
C 

C  STARTING    LAT    AND    LONG     FOR    SEARCH 

C  LAT    =    37/05/00.0  LONG    =    122/05/00.0 


C 


SLAT    =    133200,0 
SLGN    =   439500.0 


c 

c 

LATITUDE 

c 

J    =    0 

DO    10    IC    =    1,26,1 

SLAT    =    SLAT    -    6  0.0 

RULAT    =    SLAT    +       30.0 

RLLAT    =    SLAT    -       30.0 

10 

IF  (HLAT    .LT.     RULAT    -AND 

15 

J    =    IC 

C 

C 

LONGITUDE 

HLAT    .GE.     RLLAT)     GO    TO    15 


1    =    0 

DO    20    IC   =    1,16,1 
SLON    =    SLON    -    60.0 
RLLON    =    SLON    +      30.0 
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RRLON    =    SLON    -       30.0 
20  IF(HLON    .LT.     RLLON    .AND.    HLON    .GE.     HRLON)     GO    TO    25 

25  I   =   IC 

RETURN 

END 

C== ======================== ================================: 

SUBROUTINE    INV2R1 ( A, F INV, GLAT1 , GLON 1 ,GLAT2,GLON2 ,FAZ, 

1    BAZ,S) 
C 

See    Appendix    B    for  subrcu-ine   INVER1. 
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DATA    SET    EXAMPLE    -    PROGRAM    ASFSEL 


TABLE    FOR    9940-Y 


50» 


1220    /    00* 


US* 


0.56    -0.51 
-0.39 


-0.41 

-0.53  -0.44 

-0.46    -0.40 

-0.54    -0.39    -0.55    -0.45 

-0. 63    -0.57    -0.49 

-0.54    -0.50 

-0.62 


40' 


35' 

360 
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TABLE    FOR    9940-Y 
1210    /   55* 


1210    /    50' 


50" 


-0.45    -G.U7  -0.41    -0.37 

45*     -0.51    -0.50 

-0.38    -0.49  -0.50    -0.46 

-0.54    -0.51  -0.49 

-0.54    -0.57  -0.49 

-0.55    -0.55  -0.50 

40*     -0.67    -0.61  -0.50 

-0.64  -0.65 


-0.60 

-0.43    -0. 44    -0.47    -0.44 

-0.52    -0.46 


35' 
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